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I.  SUMMARY  AND  CONCLUSIONS 


1.1  THE  EXPERIMENT 

During  the  summer  of  1964  a  SOFAR  velocity  measurement  experiment  was 
conducted  by  GM  Defense  Research  Laboratories  between  Hawaii  and  *he  California 
coast  and  between  Hawaii  and  certain  islands  in  the  North  Central  Pacific.  Signals 
also  were  received  and  utilized  for  velocity  computation  from  an  independent, 
concurrent  experiment  in  the  Aleutians.  GM  DRL  developed  and  used  a  depioyabie 
acoustic  monitoring  buoy  station  to  record  signals.  Additional  fixed  acoustic 
Installations  were  used  on  the  coast  of  California  and  on  the  islands  of  Eniwetok 
and  Midway.  The  charges  for  the  GM  DRL  Oahu  shots  were  M-3  demolition 
blocks,  equivalent  to  a  few  pounds  of  TNT.  These  were  electrically  fired  at  the 
approximate  depth  of  the  SOFAR  axis  from  USS  ARKARA  (ATF  98).  In  the 
Aleutians,  shots  equivalent  to  about  one  ton  of  TNT  were  fired  at  shallow  depth 
by  the  RV  SEA  SCOPE,  and  248-lb  depth  charges  to  be  located  acoustically  were 
dropped  by  a  Navy  patrol  aircraft. 

1.2  THE  OBJECTIVES 

The  principal  objective  of  the  experiment  was  to  measure  the  average  SOFAR 
velocity  between  widely  spaced  termini  and  to  determine  the  fluctuations  in  this 
measurement.  A  secondary  objective  was  to  determine  whether,  over  a  period 
of  five  days,  coherent  temporal  variations  could  be  observed  and,  if  so,  to 
correlate  them  with  possible  causes.  A  knowledge  of  the  SOFAR  velocity  and  its 
fluctuations  is  required  for  acoustical  source-location  sysbms.  (T  phase  sources 
are  now  being  located  acoustically,  and  potential  systems  to  locate  the  sources 
of  nuclear  events  by  means  of  their  underwater  acoustic  signal*  are  the  subject 
of  current  research).  In  a  program  for  computing  acoustical  source  locations, 
the  mean  acoustic  propagation  velocity  is  an  important  parameter;  also  important 
are  the  fluctuations  of  the  acoustic  propagation  velocity.  These  fluctuations,  which 
are  a  measure  of  the  uncertainty  in  our  knowledge  of  the  propagation  velocity, 
largely  determine  the  obtainable  precision  in  the  location  of  source  events. 
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1.3  RESULTS 

The  average  propagation  velocity,  computed  and  observed  for  various  paths, 
together  with  fluctuations,  are  indicated  below: 


Source 

Location 

Receiver 

Location 

Propagation  Velocity  (meters  per  second) 

Computed 

Number  of 
Observations 

Observed 

Standard 

Deviation 

Oahu 

San  Miguel  Is. 
Calif,  coast 
(buoy  station) 

1480.87 

50 

1479.27 

0.85 

Aleutian 

Trench 

Same 

1472.87 

7 

1468.55 

0.51 

Oahu 

Midway 

1482. 09 

22 

1480.66 

0.31 

Oahu 

Eniwetok 

1472.87 

20 

1482.69 

0.51 

The  predominant  part  of  the  standard  deviations  quoted  above,  particularly  in  the 
first  two  paths,  must  be  attributed  to  experimental  effects  and  not  to  the  fluctuations 
of  the  properties  of  the  medium.  It  was  not  possible  to  correlate  the  observed 
fluctuations  in  average  propagation  velocity  with  any  geophysical  effects;  nor  was 
it  possible  over  a  five  day  period  to  determine  unambiguously  any  trends  in  the 
average  propagation  velocity.  Hence,  these  fluctuations  must  be  considered  random 
and  attributed  to  the  statistical  nature  of  the  generation,  propagation,  and  measure¬ 
ment  process. 

A  source  location  program  was  developed  and  applied  in  a  real  problem  of  locating 
the  sources  of  signals  from  the  248-lb,  air-dropped  Aleutian  charges.  An  rms 
ensemble  residual  of  about  4  km  was  observed;  the  mean  contradiction  between 
GM  DRL  acoustic  positioning  and  the  navigational  information  from  the  aircraft 
\yas  11.9  km. 

1.4  CONCLUSIONS 

If  the  standard  deviations  reported  above  are  combined  with  other  data,  the 
minimum  uncertainty  in  range  determination  may  be  calculated  for  various 
propagation  paths: 
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Paths 

Effective 
Std.  Dev. 
(m/sec) 

a 

<c> 

Average 

Range 

(km) 

Range 

Uncertainty 

(km) 

Oahu  to  West  Coast 

0.775 

5  238x10' 4 

-3884 

2.03 

Aleutians  to  West  Coast 

0.955 

6.  495  x  10'4 

4295 

2.79 

Oahu  to  M  idway 

0.310 

2.  094  x  10'4 

2151 

0.45 

Oahu  to  Eniwetok 

0.350 

2.361x10"4 

4344 

0.47 
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II.  INTRODUCTION 


2. 1  STATEMENT  OF  PROBLEM 

The  detection,  classification,  and  location  of  underground  or  underwater  nuclear 
explosions  must  utilize  a  portion  of  the  energy  of  the  detonation,  generally  after 
the  energy  has  been  propagated  long  distances  (from  the  source  of  the  blast  to 
the  location  of  the  detection  devices).  Most  systems  for  these  purposes  are 
seismic,  i.  e. ,  they  measure  energy  that  has  traveled  through  the  earth's  litho¬ 
logical  material.  Recently,  however,  there  has  been  an  increase  in  attention  to 
underwater  acoustical  methods  of  detection.  For  land-based  blasts  these  methods 

depend  on  the  T  phase,  a  low-frequency,  underwater  acoustic  signal  seismicaily 
Cl  2  3V* 

generated.  ’  ’  ’  For  in-water  blasts  the  acoustical  energy  is  generated  directly 
by  the  blast  mechanism  and  propagated  long  distances.  One  distinctive  characteristic 
of  the  in-water  blast  is  that  it  sometimes  produces  a  pronounced,  closely- spaced 
line  structure  that  is  not  seen  in  T  phases.  Source  localization  depends  on  the 
stability  of  the  ocean  as  a  long-range  propagation  medium.  The  two  principal 
objectives  of  this  experiment  were  (1)  to  obtain  a  measurement  of  the  average 
propagation  velocity  in  the  Northeast  Pacific  over  various  long  path  lengths 
based  on  a  statistically  significant  number  of  readings  and  (2)  a  measurement  of 
the  stability  of  this  velocity  over  a  period  of  several  days, 

2.2  STATEMENT  OF  V/QRK 

The  contract  statement  of  work  provided  for  the  following: 

”  (5)  conduct  field  experiments  related  to  the  generation,  propaga  ion, 
and  detection  of  T  phase  signals, 

(6)  at  designated  stations  the  Contractor  will  monitor  the  1964 
VELA  UNIFORM  Aleutian  Shot  Series . , 

RaisetT  nuinEers  in  parentheses  refer  to  list  of  references  on  page  51. 
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(7)  provide  technical  personnel  at  the  stations  for  the  purpose  of  making 
the  recordings;  and 

(8)  reduce  the  data  to  tubular  form;  and 

(9)  perform  special  analysis  on  selected  signals;  and 

(10)  provide  a  copy  or  the  original  on  magnetic  tape  of  the  recordings 
taken. " 

2.3  CONCEPT  OF  THE  EXPERIMENT 

The  SOFAR  calibration  experiment  was  an  attempt  to  measure  properties  of  the 
ocean.  Shots  were  to  be  detonate  i  from  a  firing  vessel  near  Oahu  and  received 
at  stations  along  the  California  coast  and  in  the  North  Central  Pacific.  To  remove 
the  masking  effect  of  experimental  error  it  was  necessary  to  impose  very  precise 
control  on  both  the  generation  and  i  ;ception  of  the  acoustic  signals  -  that  is,  a 
precise  knowledge  of  both  the  time  and  location  of  the  sound  sources  was  demanded. 
Since  a  statistically  significant  sample  size  was  sought,  about  100  shots  we^e 
planned.  Economy  and  the  required  precise  timing  of  the  shots  prompted  the  use 
of  Navy  M-3  demolition  blocks  (Composition  C),  electrically  detonated  ar.d 
timed.  For  the  firing  system  it  was  necessary  to  design  and  construct  an 
electrical  firing  box  with  special  timing  terminals  and  to  procure  special  seismic 
blasting  caps,  pressure  rated  for  1200  psi.  The  requirement  for  precise  knowledge 
of  the  shot  position  meant  that  the  firing  vessel  had  to  have  a  sophisticated  navigation 
control  system  -  this  was  Hydrodist,  an  electronic  distance-measuring  system, 
comprised  of  two  shipboard  units  and  two  shore  stations  located  accurately  with 
respect  to  triangulation  points. 

A  self-recording  acoustic  buoy  system  was  developed  and  used  in  the  experiment. 
This  system,  cor  listing  of  a  submerged  shallow  buoy  suspending  a  deep  hydro¬ 
phone  on  its  mooring  line  at  the  SOFAR  axis,  recorded  signals  for  about  two 
weeks.  The  system  served  two  purposes:  (1)  it  backed  up  and  supplemented  data 
from  the  fixed  stations,  and  (2)  it  tested  the  concept  of  a  deployable  recording 
buoy  which  could  be  implanteo  in  strategically  important  but  inhospitable 
geographic  areas. 
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2.4  CHRONOLOGY 

The  significant  events  in  the  conduct  of  this  experiment  are  listed  chronologically 
below  (all  dates  1964). 

4/20  Commenced  ordering  hardware  for  field  experiment 
5/18  Commenced  work  as  per  task  statement 
6,  22  Conducted  site  reconnaissance 

6/26  Liaison  established  with:  COMSERVPAC,  CO  ARKARA, 

US  ARMY  HAWAII,  University  of  Hawaii,  NAD  OAHU, 

Eaves  &  Meredith 

7/1  Design  firm  on  buoy  and  explosive  source  systems 

7/8  Pressure  tested  buoy  (MS  15-64) 

7/14  Arrangements  completed  with  OU/CWSF  for  manning  site 

7/20  Conducted  .’st  major  system  sea  tests  (MS  16-64) 

7/22  Tested  buoy  installation,  ship  positioning  with  Hydrodist, 

explosive  rigging  and  firing  (recovery  of  buoy  unsuccessful) 

7/23  Recovered  buoy  using  support  vessel  to  assist  SWAN 

7/26  Naval  Facility  reconnaissance 

8/7 

8/10  Conducted  shallow  water  release  tests  (MS  19-64) 

8/11 

8/17  Coordinated  arrangements  with  NEL 

8/19  Group  A  NavFac  field  party  on  station 

8/21  First  NavFac  contingent  recalled  temporarily 

8/25  Hawaii  field  party  departs  for  exercise 

8/27  Deployed  buoy  system  for  recording  of  firing  runs 

8/31  ARKARA  at  sea 

Group  B  NavFac  field  party  on  station 
9/1  GM  shots  commence 

9/2  Alaskan  shots  commence 

9/5  GM  shots  complete 

ARKARA  returns  to  port 

9/6  Group  A  NavFac  field  party  returns  GM  DRL 

9/S  Buoy  recovered  (MS  26-64) 

Hawaii  field  crew  returns 
9/21  Alaskan  shots  complete 

9/21  Group  B  NavFac  field  party  returns 
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III.  ACOUSTICAL  SOURCE  SYSTEM 


3.1  GENERAL 

Sources  of  the  timed  signals  had  to  be  both  safe  and  inexpensive,  and  information 
about  the  position  and  time  of  the  source  signal  had  to  be  accurate.  Safety  was  a 
concern  because  the  decision  had  been  made  to  use  available  explosives  as  sound 
sources.  Expense  was  important  because  of  the  large  number  of  shots  to  be  fired. 

Based  upon  all  of  these  considerations,  the  sound  source  system  was  designed 

to  use  electrically-fired  standard  Navy  M-3  demolition  blocks.  The  explosive  is 

Composition  C,  a  stable  plastic  material  capable  of  taking  the  1200  psi  of  the 
(1) 

firing  depth.  'Also  capable  of  withstanding  the  high  ambient  pressure  were  the 
special  seismic  blasting  caps,  manufactured  by  DuPont  and  designated  "SSS. " 

3.2  ACCURACY  LIMITS 

The  overall  design  of  the  experiment  dictated  the  precise  knowledge  required  of 
charge  location  and  detonation  time.  Variations  in  transit  time  due  to  uncertainties 
in  the  position  of  the  charge  and  time  ol  tiring  were  to  be  kept  below  the  point 
where  they  could  mask  the  obser  vable  fluctuations  in  transit  time  caused  by 
oceanic  inhomogenieties.  A  limit  to  precision  is  set  by  the  characteristics  of  the 
received  signal  itself.  The  time  of  the  SOFAR  axis  arrival,  indicated  by  an 
abrupt  fall  in  signal  level,  in  measurable  to  about  100  ms.  If  uncertainties  in 
source  position  and  time  were  each  to  be  held  to  about  20%  of  the  signal  uncertainty, 
then  they  should  not  measurably  widen  the  statistical  spread  in  the  results.  Hence, 
source  time  uncertainty  was  to  be  held  to  20  ms,  and  position  uncertainty  to  100 
feet.  Doth  of  these  were  attainable. 

The  time  measurement  required  a  firing  circuit  that  would  record  at  the  same 
time  it  fired  the  charge.  The  position  uncertainty  requirement  led  to  precise 
electronic  equipment  to  locate  the  firing  ship.  An  electronic  distance-measurement 
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system  known  as  "Hydrodist"  was  employed;  this  navigation  system  requires  two 
accurately  determined  land-based  stations  from  which  precision  distance  measure¬ 
ments  can  be  made  to  the  ship.  With  a  specially  prepared  grid  chart,  the  geo¬ 
graphical  position  of  the  ship  can  be  determined  to  uncertainties  of  less  than  10 
feet  in  ranges  up  to  14  miles  from  the  shore  stations. 

The  charges  were  suspended  on  an  electric  cable  with  a  50- lb  weight  at  the  end 
to  reduce  wire  angle  at  the  moment  of  firing;  2, 400  feet  of  cable  was  streamed. 
Correction  of  horizontal  and  vertical  displacement  depended  on  wire  angle  at  the 
ship;  the  amount  of  change  was  arbitrarily  fixed  at  one-half  of  the  straight 
geometric  correction.  This  procedure  is  based  on  the  practice  of  oceanographers 
in  correcting  Nansen  bottle  cast  geometry: 

=  1/2  L  Sin  9  horizontal  correction 

w 

AZ  =  Lw  (1  -  1/2  Cos  8)  vertical  correction 
where  L>w  =  length  of  wire 

8  =  wire  angle  at  surface,  i.e. ,  angle 

between  vertical  and  tangent  to  wir  e. 

3.3  FIRING  CIRCUIT 

The  firing  box  was  designed  to  operate  from  the  ship's  120- vdc  system.  Safety 
interlocks  of  banana- plug  jumpers  that  were  part  of  the  firing  circuit  prevented 
the  firing  circuit  from  closing  accidentally.  (If  a  technician  keeps  an  Interlock 
with  him  while  working  on  the  firing  circuit,  the  circuit  can  be  closed  only  with 
his  knowledge  and  concurrence. )  The  firing  box  also  provided  a  voltage  tap  on 
either  side  of  a  1-ohm  resistor  to  measure  the  firing  current.  Leeds  from  this 
point  went  to  one  channel  of  a  two- channel  Sanborn  recorder  to  indicate  the 
instant  of  switch  closure;  the  other  channel  received  time  marks  from  radio 
station  WWV.  The  firing  circuit  diagram  is  shown  in  Figure  1,  and  the  firing 
box  is  illustrated  in  Figure  2. 
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CHANNEL  1  or  2 


FIRING  LEADS 
2  PIN  CANNON  FEMALE 


Figure  1  Firing  Circuit  Diagram 
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Figure  2  Firing  Box 
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3.4  SOURCE  LEVELS 

Source  levels  and  spectra,  computed  from  Weston ,  were 


where  W  =  charge  weight  in  pounds 

and  R  =  distance  from  detonation  center  in  feet 

The  time  constant,  t  is  found  to  be  0. 15  ms  when  calculated  for  a  5-lb  charge, 

0  (41 

a  10-ft  horizontal  displacement,  and  a  depth  of  2, 400  feet.'  'The  first  bubble- 

(51 

pulse  frequency  will  be  88. 6  cps.  From  Cole,  the  diameter  of  the  first  bubble 
will  be  6.  7  feet.  In  summary,  the  signal  from  these  charges  is  very  sharp, 
with  rapid  decay,  a  high  bubble  pulse  frequency,  and  a  small  bubble.  The  receiving 
electronics  was  designed  to  be  sensitive  in  the  frequency  range  of  the  first  bubble- 
pulse  frequency. 
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IV.  ACOUSTICAL  RECORDING  SYSTEM  (BUOY) 


4.1  MECHANICAL  SYSTEM 

Figure  3  shows  the  configuration  of  the  deployed  acoustical  recording  buoy 
system.  The  hydrophone  is  located  in  the  vicinity  of  the  SO  FAR  axis,  and  the 
spherical  buoy  containing  recording  instruments  and  batteries  is  located  about 
100  feet  below  the  surface.  Buoyancy  of  the  submerged  buoy  provides  tension 
to  keep  the  system  vertical  and  taut,  conditions  necessary  for  minimum  horizontal 
excursion  of  the  hydrophone.  The  clump  anchor  connects  by  means  of  a  3/8-in. , 

1  x  19  steel  wire  rope  (IWRC*)  to  the  hydrophone/release  assembly,  which  in 
turn  connects  to  the  submerged  buoy  via  a  3/8-in. ,  7-conductor,  double-helical, 
armored  well- logging  cable.  The  submerged  buoy  ties  to  a  surface  marker 
through  a  1-in.  Dacron  line  containing  three  electrical  conductors,  one  in  each 
strand.  The  release  mechanism  is  controlled  by  a  pre-aet  timer  which  discharges 
an  explosive  nut,  dropping  the  eye-bolt  in  the  lower  mooring  and  freeing  the  rest 
of  the  system  to  surface  under  the  upward  pull  of  the  submerged  buoy.  The 
anchor  and  mooring  cable  are  not  recovered.  If  the  timing  mechanism  fails,  or 
if  premature  recovery  is  desired,  the  explosive  nut  may  be  fired  from  the  surface 
via  one  of  the  electrical  conductors  in  the  Dacron  line.  The  other  two  conductors 
permit  monitoring  of  the  hydrophone. 

The  water  depth  at  the  implantnent  site  was  about  4, 200  feet,  and  the  SO  FAR 
axis  was  found  at  about  2,  400  feet.  These  dimensions  required  1, 800  feet  of 
1  x  19  mooring  cable  and  2, 300  feet  of  well-logging  cable.  The  60-inch -diameter 
submerged  buoy,  weighing  386  pounds,  produced  a  net  positive  buoyancy  of 
about  2,415  pounds  when  loaded  with  electronics  and  batteries.  The  anchor, 
consisting  of  nine  discs  and  a  nose  cone,  weighed  about  5, 300  pounds  in  water. 
Since  the  submerged  buoy  supported  a  suspended  weight  of  1, 210  pounds,  a  net 
tension  of  1,  205  pounds  remained  in  the  mooring  at  the  anchor. 

*  Independent  Wound  Posin  Core 


12 


OM  DIM  Nil  IIIIANCN  LAHQIATOIISI  $>  IINIRAi.  MOTOM  CORPORATION 


TR85-28 


SURFACE 


Figure  3  Acoustic  Recording  buoy  System  Deployed 
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The  hydrophone,  a  BaTiOg  ceramic  cylinder  resonant  at  28  kc  with  a  calculated 
sensitivity  of  about  -96  dbv  per  /ibar  and  a  capacitance  of  0.0182  Mf,  is  shown 
in  Figure  4.  The  hydrophone/release  assembly  housing  is  shown  in  Figure  5. 

The  spherical  buoy,  in  position  to  be  laid  from  the  afterdeck  of  SWAN,  is  shown 
in  Figure  6. 

4.2  ELECTRONIC  SYSTEM 

Figure  7  is  a  block  diagram  of  the  electronic  system.  The  hydrophone  drives 
a  40-db  preamplifier  which  in  turn  drives  the  cable  to  the  submerged  instrument 
buoy.  The  frequency  band  at  the  output  of  the  preamplifier  is  designed  to  lie 
between  50  cps  and  500  cps,  the  lower  cutoff  being  above  the  high-level  natural 
T  phase  noise,  and  the  upper  cutoff  above  expected  signal  frequencies. 

At  the  submerged  buoy  the  signal  from  the  preamplifier  is  amplified  further, 
and  the  detected  output  drives  an  optical  galvanometer  which  marks  light-sensitive 
paper  carried  on  the  drum  of  a  helical  recorder.  The  recorder  advances  at  the 
rate  of  one  turn  every  15  minutes;  at  the  end  of  12  hours,  on  command  from  a 
crystal-controlled  timer,  the  drum  re-sets  itself  and  advances  its  chart  paper. 
Enough  paper  is  carried  for  15  days  of  continuous  recording.  The  crystal- 
controlled  timer  also  drives  a  separate  optical  galvanometer  which  places  minute 
and  hour  marks  on  the  record. 

Figures  8  and  9  show  circuit  schematics  of  the  hydrophone  preamplifier  and 
the  buoy  amplifier. 

4.3  INSTALLATION  AND  RECOVERY 

Implantment  of  this  single- taut- line  buoy  system  utilizes  the  free-fall,  baled- 
anchor  technique  developed  at  GM  DRL.  This  procedure  permits  installation  in 
a  precise  geographic  spot,  precludes  the  lengthy  handling  of  heavy,  pendulous 
weights,  and  may  be  accomplished  smartly  in  very  rough  seas. 
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RESONANT  FREQUENCY  28  KC 
CALCULATED  SENSITIVITY -96  dbv/fib 
CAPACITANCE  0.0182//{ 

PRESSURE  CAPABILITY  >2000  psi 


Figure  4  Hydrophone  Assembly 
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Figure  6  Submersible  Instrument  Buoy  Aboard  R/V  SWAN 
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Figure  7  Electronic  System  Block  Diagram 


Figure  8  Hydrophone  Preamplifier  Circuit  Schematic 
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Figure  9  Recording  Electronic  Circuit  Schematic 
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Figure  10  illustrates  the  Installation  sequence.  First  into  the  sea  goes  the  safety 
float  which  prevents  loss  of  the  system  in  case  of  trouble  during  implantment; 
the  float  is  recovered  after  a  satisfactory  installation.  Next  to  be  streamed  is 
the  Dacron  surface  pennant  which  is  supported  along  its  length  by  plastic  pneu¬ 
matic  floats.  Then  the  spherical  instrument  buoy  goes  over.  At  this  point  the 
buoy  is  floating;  later  it  will  be  carried  down  by  the  anchor.  Then  with  the  ship 
proceeding  at  about  one  knot  to  keep  the  buoys  separated  and  tension  in  the  lines, 
the  well-logging  cable  Is  streamed  until  it  is  payed  out.  The  end  of  the  cable 
attaches  to  the  upper  end  of  the  hydrophone/release  assembly,  which  is  lashed 
to  the  transom  and  ready  to  be  released.  One  end  of  the  1  x  19  mooring  cable 
attaches  to  the  underside  of  the  hydrophone/release  assembly,  and  the  other  to 
the  bale  on  the  anchor  assembly,  which  is  rigged  outboard  ready  to  drop.  The 
mooring  line  is  pre-cut  to  the  Kr.own  water  depth  at  the  installation  site.  Navigation 
of  the  ship  at  the  point  of  drop  must  be  very  precise,  especially  if  the  bottom  is 
uneven  or  steeply  sloping.  During  this  installation  by  the  SWAN,  the  Hydrodist 
electronic  survey  instruments  inboard  ship  were  used  in  conjunction  with  associated 
instruments  at  shore  stations  of  known  position.  When  the  ship  is  over  the  precise 
geographical  point  and  the  sight  in  the  instrument  cable  is  correct,  the  anchor  is 
dropped.  It  proceeds  straight  down,  paying  out  wire  as  it  goes.  The  hydrophone/ 
release  assembly  is  then  cast  overboard,  thus  freeing  the  entire  system  from  the 
Installing  ship.  When  the  mooring  wire  is  out  of  the  bale,  the  spherical  buoy 
begins  to  move  toward  the  drop  point  and  then  abruptly  submerges.  Finally,  the 
anchor  strikes  bottom  and  the  system  is  configured  as  shown  in  Figure  3,  except 
for  the  safety  float  which  is  cast  free  and  recovered.  After  the  hydrophone  is 
monitored  for  proper  output,  the  surface  pennant  is  taken  up  and  secured  to  a 
plastic  marker  float.  The  taut-line  acoustic  monitoring  system  is  now  deployed 
and  operating. 

The  helical  recorder  and  timer  are  started  before  the  installation,  when  initial 
time  calibrations  are  inserted  on  the  record  using  signals  from  radio  station  WWV. 
After  recovery,  a  terminal  time  calibration  will  also  be  recorded.  These  time 
marks  permit  any  long-term  drift  of  the  crystal-controlled  oscillator  to  be 
measured  so  that  timing  corrections  may  be  made  in  the  data. 
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Normal  recovery  is  effected  by  discharging  the  explosive  nut  mentioned  in 
Section  4. 1  through  the  circuitry  described  in  Section  4.  2.  The  actuating  signal 
may  come  either  from  a  pre-set  timer  in  the  hydrophone/release  assembly  or 
from  a  manually  thrown  switch  on  the  surface.  When  the  explosive  nut  dischrages, 
the  1  x  19  mooring  cable  is  released  and  the  system  from  the  hydrophone  upward 
rises  under  the  buoyant  pull  of  the  submerged  buoy.  Since  the  anchor  and  mooring 
are  discarded,  normal  recovery  does  not  involve  the  handling  of  large  weights. 

The  two  methods  of  recovery  provide  redundancy,  but  during  the  1964  SWAN 
operation  both  methods  failed,  forcing  recourse  to  emergency  dragging.  Recovery 
was  finally  effected  by  fine  seamanship  abetted  by  good  fortune. 

4.4  RECORDS 

A  reproduction  of  a  portion  of  the  record  recovered  from  the  submerged  instrument 
buoy  is  shown  in  Figure  11.  This  record  is  taken  from  the  period  of  operation 
when  GM  DRL  SOFAR  shots  were  being  fired  hourly,  and  several  of  the  signals 
from  these  shots  are  shown.  This  was  also  a  period  of  low  background  ambient 
noise;  the  signal-to-noise  ratio  of  these  shots  lies  between  30  and  32  db.  Duration 
of  these  signals  ranges  between  9  and  10  seconds;  their  typical  sudden  termination 
times  the  energy  which  has  traveled  along  the  SOFAR  axis.  Figure  12  shows  the 
sign?!  received  from  one  of  the  TI  Aleutian  shots  (No.  22).  Since  the  recorder 
traces  the  detected  output  of  the  amplifier  in  the  buoy  within  the  dynamic  range  of 
the  electronics,  the  area  under  the  signal  curve  is  related  to  the  total  energy 
received  during  the  shot.  The  TI  shot  signal  has  a  duration  of  about  19  seconds. 

The  acoustic  monitoring  buoy  was  installed  on  27  August  1964  in  the  following 
location,  about  6  miles  south  of  San  Miguel  Island  off  the  coast  of  Southern 
California: 

Lat  33°  46'  15”  N 
Long  120°  22'  30”  W 

The  buoy  was  recovered  on  9  September  1964. 
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Figure  11  GM  DRL  SO  FAR  Signals  Received  by  Acoustic  Recording 
Buoy  System 
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V.  THE  FIELD  OPERATION 


5.1  PARTICIPATING  GROUPS 

Two  ships  and  two  field  parties  participated  directly  in  this  shot  program,  and 
one  ship  and  field  party  participated  indirectly.  A  navy  fleet  tug,  USS  ARIKARA 
(ATF  98)  (Fig.  13)  served  as  firing  vessel  for  the  GM  DRL  Oahu  SQFAR  shots. 
Providing  tactical  support  off  the  California  coast,  the  GM  Research  Vessel 
SWAN  (Fig.  14),  in  the  course  of  several  operations  at  sea,  tested  the  charge 
firing  system  described  in  Section  3,  and  tested,  installed,  and  recovered  the 
acoustic  recording  buoy  described  in  Section  4.  One  GM  DRL  field  party 
coordinated  the  Hawaiian  firing  operation,  with  one  section  aboard  ARIKARA 
and  another  operating  the  shore-based  navigation  stations.  The  other  GM  DRL 
field  party  recorded  signals  received  at  fixed  hydrophone  arrays  on  the  West 
Coast.  A  third  ship  and  field  party,  RV  SEA  SCOPE,  a  converted  navy  AM  under 
the  control  of  Texas  Instruments,  Inc.  (TI),  al though  operating  in  a  completely 
separate  and  independent  program  was  involved  indirectly  in  the  SOFAR  calibration 
experiment.  (See  Reference  6  for  a  description  of  the  TI  program.)  The  GM  DRL 
recording  field  party  detected  both  the  small  GM  Oahu  SOFAR  shots  (5  -  18  lb) 
and  the  large  TI  Aleutian  shots  (2376  and  248  lb).  Personnel  of  the  Tsunami 
Project,  Hawaii  Institute  of  Geophysics,  University  of  Hawaii,  also  provided 
data  collection  services.  They  tabulated  arrival  times  of  the  GM  DRL  Oahu  shots 
at  two  of  the  Pacific  Missile  Impact  Location  System  (MILS)  stations  and  forwarded 
data  together  with  recorded  magnetic  tape  to  GM  DRL. 

5.2  PREPARATIONS 

The  success  of  such  an  operation  requires  preparation.  Not  only  was  there  hard¬ 
ware  to  be  designed,  ordered,  assembled,  and  tested,  but  coordination  had  to  be 
established  among  the  groups  involved  and  services  provided  through  existing 
military  organisations.  The  Office  of  Naval  Research  provided  the  basic  ship 
services.  Details  of  schedule  and  logistics  were  then  worked  out  directly  with 
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Figure  13  USS  ARKARA  (ATF  98) 
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Figure  14  GM  DRL  Research  Vessel  SWAN 
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Commander  Service  Force,  Pacific  Fleet,  and  the  Commanding  Officer  of  TJSS 
ARIKARA.  To  establish  the  shore  navigation  stations,  permission  had  to  be 
obtained  through  Headquarters,  U  Army,  Hawaii.  TheM-3  demolition  blocks 
were  provided  by  the  Naval  Ammunition  Depots,  one  at  Oahu  and  the  other  at 
Seal  Beach,  Calif. ,  after  an  official  letter  request  by  the  Office  of  Naval  Research 
and  telephone  contact  by  GM  DRL.  Communication  arrangements  were  made  and 
information  exchanged  with  TI  and  Air  Force  Cambridge  Re  arch  Labs  (AFCRL), 
the  governmental  control  agency  for  the  Aleutian  shots.  Contact  was  established 
also  with  the  Navy  Electronics  Laboratory,  San  Diego,  and  their  requirements 
were  incorporated  into  the  data  recording  method. 

5.3  SCHEDULE  AND  PERFORMANCE 

The  GM  Oahu  SOFAR  shots  were  to  be  fired  once  an  hour  for  five  days,  which 
amounts  to  120  detonations.  Due  to  difficulties  of  one  kind  and  another,  only 
half  cf  this  number  were  actually  fired.  Figure  15  shows  the  actual  firing  times 
of  these  shots;  note  that  good,  around-the-clock  coverage  was  provided  with 
the  exception  oi  one  hour  when  a  series  of  misfires  thwarted  efforts. 

The  original  schedule  of  the  TI  Aleutian  shots  and  the  actual  performance  are 
shown  together  in  Figure  18. 
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Figure  15  Diurnal  Distribution  of  GM  DRL  SO  FAR  Shots 
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VI.  DATA  PROCESSING  AND  RESULTS 


6.1  GENERAL 

Appendices  I  and  II  contain  data  summaries  of  both  the  source  events  and  recorded 
signal  arrivals.  Source  data  Included  detonation  time  and  position,  charge  depth 
and  weight,  shot  number  and  qualitative  performance.  Received- signal  data 
included  arrival  times  at  two  hydrophones,  duration,  signal  level  increase, 
background  level,  and  amplifier  gain.  From  source  and  arrival  time  information 
at  a  single  receiver,  the  average  propagation  velocity  was  calculated  for  each 
shot;  from  this  Information  at  a  group  of  receivers,  Bource  locations  were 
computed;  and  from  the  information  provided  by  a  large  number  of  arrivals, 
statistics  relating  to  the  stability  of  the  SOFAR  velocity  were  derived. 

6.  2  SOFAR  VELOCITY  COMPUTATIONS 
Computer  Program 

Reduction  of  the  observed  source  and  received  tiraeB  was  accomplished  on  an 
IBM  7040  digital  computer  in  two  principal  stages.  The  firBt  corresponds  to  the 
production  of  a  source  log  deck.  The  second  used  thiB  deck,  plus  rece  /ed 
times,  to  compute  average  sound  velocity. 

3ource  Log  Program 

This  first  program  included  bringing  together  several  separate  logs  kept  by  the 
charge  detonation  field  party  and  performing  several  simple  calculations  to 
yield  a  clear  source  deck.  The  principal  stepB  were  as  follows: 

a)  Assembling  of  shot  number,  date,  time,  charge  weight,  and 
performance  from  logs 

b)  Reduction  of  ship  geodetic  position  from  recorded  Hydrodist 
readings  to  degrees,  minutes  and  nearest  second 
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c)  Correction  of  ship's  position  to  shot  position  from  logs  of  wire 
payed  out,  wire  angle,  and  true  bearing.  (In  an  attempt  to 
compensate  for  wire  catenary,  the  depth  and  horizontal  range 
of  the  shot  from  the  ship  were  calculated  using  an  effective 
wire  vertical  angle  of  one-half  the  recorded  value. ) 

d)  Printing  and  punching  of  all  these  quantities  (the  resulting 
source  log  is  presented  in  Appendix  I). 

Sound  Velocity  Program 

The  second  stage  step  included  going  from  this  source  log  deck  and  receiver  log 
decks  to  average  sound  velocity  (listings  of  the  receiver  decks  are  printed  in 
Appendix  n).  The  principal  data  given  for  each  shot  are  station  and  shot  numbers, 
date,  arrival  times  at  one  to  four  hydrophones,  shot  trace  amplitude,  and  system 
amplitude  calibration  (a  listing  of  the  sound  velocity  program  is  given  in  Appendix 
HI).  The  principal  events  are: 

a)  Input  of  alpha-numeric  constants  used  by  the  program 

b)  Reading  and  storing  of  source  deck  -  information  is  referenced 

by  shot  number  * 

c)  Conversion  of  geodetic  source  positions  to  geocentric  coordinates 
in  radians 

d)  Reading  of  geodetic  coordinates  of  centers  of  two  receiving 
hydrophone  arrays  (in  radians)  and  conversion  to  geocentric 
radians  (in  this  program  It  was  assumed  that  the  first 
coordinate  referred  to  the  received  times  on  channels  1  or 

3  on  the  Sanborn  trace,  or  both,  and  the  second  coordinate 
refers  to  channels  2  and  4). 

e)  Reading  one  receiver  card  corresponding  to  a  particular 
shot  number 
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Computation  of  great-circle  distance  in  radians  from  shot  source 
location  to  receiver  (uses  a  subroutine  entitled  GCD).  This  is 
converted  to  a  distance  in  meters  using  the  radius  of  the  earth 
at  the  mean  of  the  source  and  station  latitudes  (see  U.S.  Nautical 
Almanac  on  International  Ellipsoid  of  reference). 

Computation  of  average  velocity  from  this  range  and  the  difference 
between  source  and  received  times 

Taking  of  time  differences  in  seconds  between  the  times  from 
channels  1  and  2  or  3  and  4.  This  relates  to  the  shot’s  apparent 
bearing. 

Computation  of  approximate  peak  sound  pressure  level,  according 
to  the  iormula 

SPL  *R+B+A-G-A  +  M 
o  o 

R  is  the  rise  (in  db)  of  the  Sanborn  trace 

B  is  the  background  level  or  reference  in  db/  1  volt 

A  is  the  current  system  attenuator  setting 

Gq  is  the  system  gain  at  some  previous  time 

Aq  was  the  attenuator  setting  when  Gq  was  measured 

M  is  the  hydrophone  sensitivity  in  db/  1  volt/pbar 

All  the  above  quantities  are  read  in  for  each  :'iot,  except  M, 
which  was  assumed  to  be  -80  db/  1  volt/microbar  at  100  cps. 

The  sound  levels  thus  computed  should  not  be  relied  on  for 
absolute  value  or  for  relations  between  stations  due  to  our 
Ignorance  of  absolute  calibration.  They  are  internally  consistent, 
however. 

Outputting  of  all  calculated  data  in  printed  form  and  punched- 
card  form  for  statistical  analysis 

Repetition  of  above  steps  e  -  j  for  each  shot  at  a  station 
Repetition  of  steps  d  -  *  for  each  receiving  station 
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The  procedure  above  was  repeated  for  each  of  the  stations,  but  not  for  the  buoy 
data.  Since  only  one  hydrophone  channel  was  available  and  the  recording  medium 
different,  the  input  format  and  calculations  had  slightly  different  forms  from 
the  above. 


A  listing  of  the  computed  sound  velocities,  time  differences,  and  sound  pressure 
levels  is  given  in  Appendix  IV  for  the  GM  shots  received  at  various  stations  and 
for  some  TI  shots  received  at  fixed  installations. 


8.3  SOFAR  VELOCITY  SUMMARY 

Average  sound  velocity  results,  computed  for  the  indicated  termini  points,  are 
discussed  below. 

Source  Receiver 

GM  DRL  Oahu  Recording  Buoy  Station 

Midway  MILS  Station 
Eniwetok  MILS  Station 
TI  Aleutian  Recording  Buoy  Station 

These  time  series  data  are  plotted  in  Figures  17  thx  jugh  20.  Mean  sound 
velocities  for  each  pair  of  termini,  together  with  standard  deviations,  are  listed 
in  Table  I  which  also  shows  the  calculated  average  sound  velocity  based  on 
Johnnon.^ 


Average  expected  sou  id  velocity  over  a  total  distance  D,  where  the  theoretical 
propagation  velocity  is  known  in  N  piecewise  increments  and  assumed  to  be 
constant  over  the  distance  interval  ,  is  given  by 


<c  > 


D 

nr 

e  -t- 
i=i  s 


where 

N 

D  -E  D 

J-l  3 


(1) 


(?.) 


I 
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SOUND  VELOCITY  (meters/sec) 
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Figure  18  Average  Sound  Velocity  of  GM  DRL  Oahu  SOFAit 
Shots  Received  by  Midway  MILS  Station 


X 

i 


37 


f  *  OKFKNSK  RiSiAitCN  LAtOMTODIIl  9  OtNKR 


*L  motors  CORPORATION 


t 


SM  OCFSNS&  KIIEUCH  U  lOK  &TOAI  IS  $  S  t  '!  A  M  MOTORS  CORPORATION 


TR65-28 


If  Uie  total  distance  D  is  divided  into  N  increments  of  equal  length,  ,  in 
which  the  SOFAJR  velocity  is  approximated  by  C.  ,  Equation  (1)  reduces  to 


<c  >  — 


N 

NL 


i=l 


-1 


(3) 


The  measured  velocity  from  Oahu  to  the  buoy  and  other  west  coast  receivers  is 
slightly  less  than  the  computed;  the  decrement  is  about  1  meter  per  second. 

(71 

Possible  reasons  for  this  effect  are  many.  In  the  first  place,  Johnson's  chart,  ' 
which  served  as  the  basis  of  the  computation,  was  in  rather  small  scale  for 
these  purposes;  hence,  the  uncertainty  of  the  computed  velocity  may  be  1/2 
meter  per  second,  or  so.  Another  hypothesis  assumes  that  the  sound  actually 
traveled  farther  than  the  strict  geodesic  calculation  would  indicate.  Variations 
in  path  could  be  in  the  horizontal  plane,  due  to  horizontal  velocity  gradients  or 
scattering,  or  in  the  vertical  plane,  due  to  changes  in  depth  of  the  SOFAR  axis. 
Such  an  effect  might  occur  if  significant  internal  wave  activity  were  present.  To 
check  the  horizontal  deviation  hypothesis,  statistical  correlations  were  run  at 
each  West  Coast  installation  between  the  sound  velocity  and  azimuth  of  arrival. 
Results  were  negative,  i.  e. ,  there  was  no  significant  correlation  between 
deviations  from  the  great  circle,  as  received,  and  the  average  sound  velocity. 
The  effect  of  vertical  variation  in  SOFAR  axis  depth  seems  to  be  in  the  right 
order  of  magnitude  and  could  account  for  reductions  in  apparent  sound  velocity 
of  the  magnitudes  observed.  The  table  below  indicates  the  relationship  between 
wavelength  and  amplitude  to  produce  a  velocity  decrement  of  one  meter  per 
second  and  ten  meters  per  second,  assuming  a  sinusoidal  path. 


Table  I 

PATH  AMPLITUDE  AND  WAVELENGTH  PRODUCING 
A  GIVEN  VELOCITY  DECREMENT 


Velocity  Decrement 

Velocity  Decrement 

1  ip/sec 

10  m/sec 

Amplitude 

Wavelength 

Wavelength 

(meters) 

(meters) 

(meters) 

10 

1200 

380 

25 

3000 

950 

50 

6000 

1900 

_ 1 
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Based  on  the  indicated  results  at  the  buoy  station  and  others,  the  net  standard 
deviation  at  a  range  of  about  2115  nautical  miles  (3923  km)  is  about  0.  05  km 
per  second.  This  uncertainty  in  velocity  Droduces  a  position  uncertainty  in  the 
source  of  about  2.  25  kilometers.  This  statement  of  uncertainty  is  one  of  the 
main  objectives  of  the  experiment.  The  remainder  of  the  receivers,  particularly 
the  MILS  station,  showed  a  smaller  spread  in  the  data,  which  would  indicate  a 
smaller  position  u  ertainty. 

6.4  SOURCE  LOCATION 
General 

During  the  course  of  the  Aleutian  exp^iment  several  248-lb  bombs  were  dropped 
by  Navy  P2V  patrol  planes. 

A  ’.urate  source  locations  were  required  for  the  TI  seismic  experiment;  accord- 
ingly  ARPA  requested  GM  DRLV  ;to  locate  the  source  of  these  shots  acoustically 
as  a  check  and  as  a  supplement  to  the  navigational  information  provided  by  the 
aircraft.  A  similar  request  was  sent  to  the  Hawaii  Institute  of  Geophysics, 
University  of  Hawaii.  The  receiving  stations  in  the  GM  DRL  source  location 
computational  net  were  spaced  along  the  west  coast  of  the  United  States  .\nd 
throughout  the  North  Central  Pacific.  Approximately  six  to  eight  station  times 
were  involved  in  each  location  computation.  Stations  used  by  the  University  of 
Hawaii  included  those  in  the  North  Central  Pacific  but  only  one  on  the  west 
coast  of  the  United  States  mainland. 

Scurce  Location  Program 

Given  the  arrival  time  of  signals  from  an  acoustic  event  at  three  or  more  known 
locations  on  the  earth's  surface,  source  time  and  location  are  to  be  determined. 
Because  signal  waveforms  become  distorted  during  propagation  it  is  difficult 
to  measure  arrival  times  precisely,  and  errors  will  appear  in  the  measured 
data.  Furthermore,  the  average  propagation  velocity  between  source  and 
receiver  may  not  be  known  exactly.  Thus,  for  more  than  three  receiving  stations 
it  is  quite  likeL-  that  there  is  no  source  location  on  the  earth's  surface  from 
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which  the  computed  and  observed  arrivals  at  all  stations  match  precisely.  This 
situation,  as  is  usual  when  attem  Hing  to  match  an  analytic  expression  to  the 
real  world,  required  curve  fitting  with  a  least  squares  criterion.  The  source 
location  and  time  are  chosen  to  minimize  the  rms  deviation  of  the  set  measured 
arrival  times  from  those  predicted  from  the  solution. 

Because  of  the  non-linearity  of  the  spherical  geometry  no  exact  analytical 
solution  is  possible  for  the  minimization  problem.  Good  results  may  be  obtained, 
however,  by  linearizing  the  equations  through  a  perturbation  approach,  i.e. , 
solving  for  small  corrections  to  approximate  solutions.  With  sufficient  iteration 
any  degree  of  precision  may  be  attained  in  finding  the  optimum  point  consistent 
with  the  data.  Note  that  bad  input  data  may  cause  this  point  to  be  an  inaccurate 
solution.  The  only  claim  is  that  the  position  is  the  best  that  the  data  supports. 

In  the  following  analysis  a  good  first  guess  for  source  time  and  location  will 
be  assumed.  The  actual  limitations  on  this  first  approximation  for  convergence 
will  be  discussed  later. 


Basic  Equations 

The  j*  station  residual  is  defined  as  the  difference  between  the  observed  arrival 

time  Tn.  ,  and  the  computed  arrival  time  T  .  . 

Uj  cj 


s. 

J 


T  .  (  L 
cj 


Tn) 


(4) 


where 

L  =  longitude 
X  =  latitude 


The  computcu  arrival  time  is  a  function  of  the  source  location  and  time.  For  a 
set  of  n  stations  the  ensemble  residual  is  defined  as  tire  rms  of  the  individual 
station  residuals. 


S  = 


(5) 
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The  source  location  problem  Is  to  find  the  geographical  point  and  origin  time 
that  will  minimize  the  ensemble  residual,  This  is  effected  through  maximizing 
the  negative  change  in  S  from  the  trial  point,  assuming  a  planar  geometry. 
The  change  in  S  may  be  expressed: 


Maximizing  the  series  in  parentheses  will  find  an  approximate  solution  to  the 
optimum  AS. 

In  Figure  21  the  initial  trial  point  is  (L  ,  X  ,  Tq  )  and  the  new  point  is 
(L  ,  X  ,  Tq)  .  Plane  geometry  and  long  distances  to  receiving  stations  are 
assumed.  With  this  assumption,  we  set 


AY  =  al 
AX  =  Cos  L  AX 


(7) 
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It  can  be  shown: 


ST  . 

Tl  4L 


Ay  Cos  « 
- I 


3T  .  AX  Sin  fi. 


(8) 

(9) 


3T  I 

wjUT<r  *to 


(10) 


The  Job  is  thus  to  find  the  AX  ,  AY  ,  and  ATq  which  will  minimize  ^  ■  ASjj' 

with  the  set  of  equations 


AX  Sin  +  AY  Cos  ^  +  C  ATQ  =  C  ASj 

AX  Sin  fi  +  AY  Cos  /S  +  C  A  TQ  =  C  As^  (11) 

A X  Sin  3  +  AY  Cos  fi  +CATn  =  CAs 
n  nun 


If  there  are  but  three  receivers,  there  are  three  equations  in  the  set  (11),  and 
an  exact  solution  to  s^  =  0  is  possible,  within  certain  geometric  constraints. 

For  four  or  more  receivers  a  least  squares  criterion  is  required.  After  several 
iterations  AX,  AY,  and  ATq  approach  zero  as  S  approaches  its  minimum. 

A  solution  is  considered  reached  when  AX  and  AY  are  both  less  than  10  meters. 


Program  Mechanics 

Appendix  III  contains  a  listing  of  the  GM  DRL  source  location  program  specially 
configured  to  solve  for  the  TI  248-lb  shot  locations  using  data  largely  supplied 
by  and  in  the  format  of  the  University  of  Hawaii  epicenter  location  program.  The 
principal  stepe  in  the  program  are  as  follows: 

a)  Alpha-numeric  constants  are  read  in. 

b)  Hydrophone  station  location  data  (in  geodetic  radians)  and  sound 
velocities  in  m/sec  appropriate  to  each  station  arc  read  in  and  filed  by  number. 
Latitudes  are  converted  to  geocentric  radians.  (Average  sound  velocities  to  the 
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West  Coast  stations  were  measured  from  the  TI  2376-lb  shots.  MILS  station 
velocities  were  supplied  by  the  University  of  Hawaii. ) 

c)  A  common  starting  point  for  all  iteration  is  read  in  and  converted 
to  radians. 

d)  A  reading  of  the  time  of  day  is  taken  from  the  real-time 
computer  clock. 

e)  Arrival  times  of  one  shot  at  each  station  that  received  it  are 

read  In,  and  the  corresponding  station  locations  are  taken  from  general  storage 
and  placed  in  working  arrays.  Arrival  times  are  converted  to  seconds. 

f)  The  variables  L  and  X  are  given  the  initial  values  in  radians  read 
in  at  c)  above.  The  variable  T^  is  given  a  starting  value  1688  seconds, 
earlier  than  the  first  arrival  time  read  in. 

g)  The  great-circle  ranges  from  the  trial  point  to  each  station  are 
computed  using  a  subroutine  GCD.  The  subroutine  also  returns  the  values  of 
the  sine  and  cosine  of  the  true  bearing  of  each  station  from  the  trial  point. 

h)  A  sound  velocity  in  radians/sec  is  computed  from  the  station  sound 
velocity  and  the  average  radius  of  the  earth  over  the  path. 

i)  The  station  residuals  are  formed,  Equation  (4). 

j)  A  subroutine  LSTSQ3  is  called.  This  solves  the  system  of 
Equations  (11)  using  the  coefficients  computed  in  g)  and  h)  for  the  values  of 
AX  ,  AY  ,  and  A T  to  minimize  the  sum  of  the  squares  of  the  residuals.  The 
subroutine  also  returns  the  resulting  value  of  the  rms  residual,  Equation  (5). 

k)  The  coordinates  L,  X,  and  T^  are  corrected  with  AX  ,  Cos  L  , 
AY,  and  At  ;  X  and  L  are  reduced  modulo  2  it  and  to  a  reasonable  latitude 
respectively. 

l)  If  the  correction 3  AX  and  AY  are  greater  than  10  meters  thi 
program  repeats  steps  g)  -  k)  until  they  are  both  10  meters  or  less. 

m)  When  the  solution  is  finished,  the  locations  are  reduced  to  geodetic 
degrees,  minutes,  and  seconds,  and  the  time  to  hours,  minutes,  and  seconds. 
Residuals  are  converted  to  meters.  A  heading,  the  final  values  of  all  residuals, 
and  the  source  coordinates  are  printed  out.  The  source  coordinates  are  also 
punched  out  on  cards. 
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n)  A  second  clock  reading  is  taken  and  the  total  elapsed  time  is 
computed  and  printed  out. 

o)  Steps  d)  -  n)  are  repeated  for  ail  shots. 

In  the  production  runs  for  the  TI  248-lb  shots,  the  above  calculations  were 
run  for  6-13  hydrophones  at  4  -  5  stations  to  a  solution  precision  of  10  meters 
in  7  -  12  iterations.  The  average  time  was  about  2.  5  seconds  for  a  complete 
solution,  including  ail  input-output  operations.  The  results  are  listed  in  Appendix 
VI,  together  with  the  results  of  similar  runs  from  the  University  of  Hawaii  (using 
MILS  data  only)  and  the  radar  fixes  on  the  plane  at  the  times  of  drop. 

Convergence  Behavior 

The  linear  set  (11)  will  always  have  a  solution;  however,  this  solution  is  not 
necessarily  the  solution  to  the  geographical  problem.  Because  the  earth  is 
spherical  and  finite,  the  assumptions  leading  to  (11)  are  never  exactly  true. 

Thus  when  the  AX  ,  AY  ,  AT  are  added  to  the  trial  parameters,  the  result 
is  still  in  geographical  error  (although  usually  an  order  of  magnitude  better), 
and  the  problem  must  be  re-soived  to  get  an  even  better  answer.  In  the  working 
version  the  Iteration  is  cut  off  when  AX  and  AY  are  both  less  than  10  meters. 
The  number  of  iterations  necessary  to  achieve  this  precision  varies  according 
to  the  accuracy  of  the  first  guess.  For  a  first  guess  within  5-10  degrees  of 
the  true  location,  5-7  iterations  are  adequate  to  reach  the  10-meters  precision. 

Allowing  a  first  guess  of  nearly  unlimited  sloppiness,  this  technique  appears  to 
have  a  very  wide  convergence  range.  The  three  probable  reasons  are  as  follows: 

1)  The  test  of  a  successful  iteration  is  whether  it  moves  the  trial 
point,  closer  or  farther  away  from  the  actual  point.  Thus,  considerable  error 
(in  fact,  ±90°)  is  still  tolerable  In  the  direction  of  movement  from  the  true  one. 

2)  The  earth's  surface  is  finite  and  closed.  Thus,  while  some  really 
wild  A  X  and  AL  can  arise  (e.g. ,  +  300°  ,  -653°),  the  trial  point  is  reduced 
modulo  360°  and  stays  on  the  surface. 

3)  The  third  variable,  time,  is  truly  linear.  Thus,  a  really  bad  first 
guess  Is  more  readily  corrected. 
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The  following  is  an  example  of  the  latitude  allowable.  Received  times  appropriate 
to  a  source  time  of  about  1800  hour  s  and  a  source  location  in  Hawaii  were  given 
to  the  program.  The  initial  guess  was  zero  hours,  zero.  Very  occasionally,  with 
a  really  bad  first  guess,  the  solution  will  converge  to  a  local  rather  than  a  global 
minimum.  The  station  residual  as  a  function  of  position  on  the  earth's  surface  is 
very  complicated  and  can  have  ore  or  more  local  minima.  As  an  example,  the 
arrival  times  of  TI  shot  74  at  the  MILS  and  West  Coast  stations  were  taken,  and 
for  a  lattice  of  points  on  the  earth's  surface  the  best  possible  source  time  was 
calculated  and  an  rms  station  residual  calculated.  Iso-residual  contours  were 
then  interpolated  and  plotted  in  (Fig.  22)  for  a  100-deg  square  in  the  Pacific  Ocean. 

As  may  be  seen,  besides  the  absolute  minimum  of  the  real  solution,  there  is  one 
other  local  minimum  near  9°N,  158°W.  Viewed  in  this  manner,  the  program 
solution  becomes  one  of  rolling  to  the  lowest  point  on  the  map.  The  weak  minimum 
did  capture  several  test  solutions.  There  is  little  danger  of  mistaking  this  solution 
for  the  true  one,  however,  since  the  residual  has  a  value  of  750  km  at  this  point. 

In  a  production  version  of  the  program  some  algorithm  should  be  used  to  leap  out 
of  such  false  minima.  One  successful  algorithm  was  to  drop  one  hydrophone  from 
the  solution  and  resolve.  This  yielded  the  correct  answer  in  every  case  tried  but 
was  expensive  in  terms  of  computer  time.  Perhaps  merely  solving  the  problem 
over  from  a  different  starting  point  would  be  best. 

The  contours  in  Figure  22  are  a  very  special  case  corresponding  to  one  particular 
receiver  configuration  and  one  particular  actual  source  location.  The  contours 
could  cnange  greatly  by  changing  the  numbers  or  locations  of  hydrophones  or  by 
changing  the  general  source  location.  A  study  consisting  of  systematic  variations 
of  these  parameters  might  prove  valuable  for  optimizing  patterns  of  receiving 
stations. 

Results  of  Computations 

Source  locations  for  the  248-lb  Aleutian  shots  by  three  agents  and  two  methods 
(GM  DRL  acoustical,  University  of  Hawaii  acoustical,  and  the  aircraft  navigational) 
ire  presented  in  Appendix  VI.  A  measure  of  the  goodness  of  fit  of  a  source 
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Figure  22  Pacific  Ocean  Area  for  Illustrative  Source  Location  Problem 
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location  calculation  to  the  data  is  the  ensemble  residual.  For  the  GM  DRL 
computation,  the  ensemble  residual  is  the  rms  result  of  the  individual  station 
residuals  expressed  in  meters.  The  station  residual  is  the  distance  the  station 
would  have  to  be  moved  toward  or  away  from  the  assumed  source  location  to 
match  the  computed  and  actual  arrival  time.  University  of  Hawaii  residuals  are 
in  seconds  of  time.  The  two  are  convertible  to  one  another  through  the  nominal 
S'Mind  velocity,  1480  meters  per  second.  Another  indicator  in  a  positional 
calculation  where  more  than  one  method  is  used  is  the  geographical  distance 
between  locations,  determined  by  separate  means.  The  averages  of  this  "distance 
contradiction"  among  the  three  location  computation  agents,  together  with  mean 
ensemble  residuals  for  the  entire  group  of  248-lb  shots,  are  Indicated  below. 

Average  Ensemble  Residuals 


GM  DRL 

4026  meters 

2. 72  seconds 

U  of  H 

3049  meters 

2.06  seconds 

Average  Location  Contradictions 

GM  DRL 

\ 

Navigation 

11.9  kilometers 

U  of  H 

Navigation 

12. 0  kilometers 

GM  DRL 

U  of  H 

7.  79  kilometers 
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Table  n 

AVERAGE  SOUND  VELOCITY 
COMPUTED  AND  OBSERVED 


Source 

Receiver 

Vp  (m/sec)  ■ 

Ve  (m/sec) 

Ve-Vp 

S 

N 

1 

3 

1480. 66 

1482.09 

1.43 

0.31 

OO 

1 

9 

-- 

1483.  32 

-- 

-- 

-- 

1 

10 

1482.69 

1484. 11 

1.42 

0.35 

20 

1 

11 

1479.27 

1480.87 

1.60 

0.85 

50 

2 

11 

1468.55 

1472.87 

4.32 

0.51 

7 

Legend. 


Receiver: 


Source: 


Vp  -  determined  from  observed  data  and  GM  DRL  program 
Ve  -  estimated  from  U  of  H  contour  plots 
S  -  standard  deviation  for  Vp 
N  -  number  of  points  used  to  determine  Vp 

8  -  Midway 

9  -  Wake  Island 

10  -  Eniwetok 

11  -  GM  DRL  acoustic  monitoring  buoy 

1  -  GM  DRL  Oahu  shots 

2  -  T1  Aleutian  Shots 
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APPENDIX  I 

SOURCE  DATA  SUMMARY 

This  appendix  includes  the  listings  of  the  final  source  logs  for  the  T1  23 ’’6 -lb 
Aleutian  shots  as  reported  in  Reference  9,  and  the  GM  DRL  Oahu  shots  used 
in  all  calculations  of  sound  velocity.  These  numbers  were  assembled  from 
various  logs  kept  at  the  time  of  the  shots.  The  reduction  is  briefly  described 
in  Section  6.  The  column  headings  are  self-explanatory. 


SOURCE  times  AND  LOCATIONS  T|  2376  Ld  SHO^S 


SHOT 

DATE 

GC 

r 

lat; tuoe 

LONGITUDE 

DEPTH 

CHAHGt. 

JfATEw 

MO  DY  YR 

HRMN 

SEC 

OG 

MN 

SC 

OLG 

MN 

SC 

METERS 

WT  Ld-* 

DEP  r  rl 

2 

23 

9/  2/64 

231  1 

00.00 

51 

29 

!  ON 

1  75 

52 

30W 

91  . 

2376. 

3->dC 

2 

22 

9/  4/64 

2336 

59.75 

51 

22 

l  5N 

1  76 

42 

22  w 

91  . 

2376. 

2104 

2 

21 

9/  b/64 

0246 

00.06 

51 

20 

l  ON 

l  77 

20 

55W 

91  . 

2  376. 

2122 

3 

20 

9/  5/64 

0354 

00.02 

51 

1  7 

43N 

l  77 

49 

55W 

91  . 

2376. 

1220 

2 

1  1 

9/  7/64 

0527 

00.00 

51 

03 

3  ON 

l  75 

51 

30  W 

91  . 

2376. 

3475 

2 

1  3 

9/  7/64 

1824 

00.02 

51 

1  4 

45N 

l  74 

19 

ocw 

91  . 

2376* 

4505 

2 

26 

9/  7/64 

2305 

00.35 

51 

45 

5  ON 

l  73 

35 

30W 

91  . 

2376. 

2415 

2 

25 

9/  8/64 

0230 

59.  75 

51 

4  1 

30N 

l  74 

14 

30W 

91  . 

2376. 

3038 

2 

24 

9/  0/64 

0557 

00.17 

51 

35 

55N 

1  75 

00 

4  5W 

91  . 

2376. 

3075 

2 

1  0 

9/  0/64 

234  1 

00.05 

51 

1  1 

Or./N 

l  78 

27 

20W 

91  . 

2376. 

24  j3 

2 

18 

9/  9/64 

0627 

00.25 

51 

02 

30N 

1  78 

55 

SOW 

91  . 

2376. 

3075 

2 

32 

9/10/64 

0522 

00.  14 

51 

4b 

42N 

1  78 

33 

4  5W 

91  . 

2376. 

1510 

2 

62 

9/1 l /64 

0931 

59.92 

52 

i  2 

00N 

l  78 

05 

40W 

61  . 

2376. 

3185 

2 

45 

9/12/64 

0155 

59.89 

52 

08 

OON 

1  79 

01 

20W 

61  . 

2  376. 

2965 

2 

42 

9/ .2/64 

0345 

00.00 

51 

51 

1  5N 

1  79 

33 

oaw 

61  . 

2  376. 

1757 

2 

40 

9/12/64 

2030 

59.92 

51 

51 

50N 

1  78 

34 

56E 

61  . 

i  376. 

3d6 

2 

41 

9/13/64 

0201 

00.15 

51 

40 

15N 

1  79 

28 

1  bw 

61  . 

2376. 

1244 

2 

30 

9/13/64 

0605 

59.97 

51 

21 

55N 

1  79 

47 

3  5W 

61  . 

2376. 

1538 

2 

46 

9/15/64 

1906 

00.  20 

52 

20 

22N 

l  77 

16 

2  3W 

91  . 

£376. 

3335 

2 

47 

9/17/64 

2206 

CO.  1  3 

52 

21 

30N 

l  76 

42 

35  W 

91  . 

2376. 

34  75 

2 

59 

9/18/64 

0122 

00.21 

52 

4b 

1  3N 

l  76 

4  3 

30W 

9  i  . 

2376. 

3550 

2 

54 

9/18/64 

1904 

00.02 

51 

28 

20N 

1  79 

23 

1  OW 

91  • 

2376. 

1061 

2 

1  7 

9/18/64 

2256 

00.23 

51 

02 

!  2N 

179 

44 

SOW 

91  . 

2376. 

3  755 

( 


1-1 


SM  DEFENSE  RESEARCH  LABORATORIES  $  GENERAL  MOTORS  CORPORATION 


TR65-28 


SOURCE  TIMES  AND  LOCATIONS  GM/DPL  SHOTS 


SHOT 

da 

GCT 

L A  T 1  TUUL 

LONG  1 TUDt 

DEPTH 

CHAHLL 

lxplosicn 

MO  DY  Y9 

HHMN 

SEC 

DG 

MN 

SC 

DEG 

MN 

sc 

METtWS 

WT 

L8j 

st-jengtm 

1 

1 

9/ 

1  /t>4 

Obo  l 

Ob 

00 

21 

4b 

2aN 

168 

1  0 

1  /  w 

72  0. 

4 

30 

GOOD 

1 

2 

i / 

1  /64 

1  000 

32 

9  0 

21 

4b 

0  ON 

1  Db 

10 

30  W 

720. 

4 

30 

GUoO 

1 

b 

9/ 

1  /64 

1600 

6  b 

60 

2  1 

17 

55N 

1  3b 

1  3 

2bW 

649. 

4 

30 

GOOD 

1 

0 

9/ 

1  /64 

1U00 

4  / 

40 

21 

4  7 

1  9N 

I  bb 

1  1 

9W 

732  . 

4 

30 

GOOD 

1 

7 

)/ 

1  /64 

200  3 

20 

00 

21 

47 

39N 

1  58 

1  0 

4dW 

70  7. 

4 

30 

GOOD 

1 

b 

9/ 

1  /64 

2201 

00 

70 

21 

47 

2  7N 

1  38 

1  0 

52  W 

676. 

4 

30 

GOOD 

1 

9 

9/ 

2/64 

0002 

64 

80 

21 

47 

26N 

1  38 

1  1 

4W 

94  7, 

4 

50 

GOOD 

1 

in 

9/ 

2/64 

"'boo 

00 

07 

21 

40 

2N 

1  68 

1  0 

42W 

711. 

4 

30 

GOOD 

1 

l  l 

9/ 

2/64 

0500 

32 

00 

21 

4H 

8N 

1  38 

10 

55W 

879, 

4 

30 

DUD 

t 

12 

9/ 

2/64 

1  100 

53 

70 

21 

47 

56  N 

1  68 

1 1 

7  W 

8  72. 

4 

30 

GOOD 

1 

1  J 

)/ 

2/64 

1  JO  3 

2b 

20 

21 

48 

IN 

1  58 

1 1 

3  W 

801  . 

4 

30 

GOOD 

1 

1  -3 

9/ 

i/64 

1  700 

35 

00 

21 

48 

1  IN 

1  58 

10 

4  1  W 

771  . 

4 

30 

GOOD 

1 

1  6 

9/ 

2/64 

1  9  0  0 

30 

00 

21 

48 

2N 

1  5  b 

1  0 

46W 

744  , 

4 

oO 

CAPS 

1 

1  / 

'-/ 

2/64 

2  1  00 

2  7 

70 

21 

47 

4  ON 

1  58 

1 1 

7W 

676. 

4 

30 

CAPS 

1 

Id 

9/ 

2/64 

230  2 

34 

/a 

21 

4b 

■jN 

158' 

10 

42W 

730. 

4 

30 

GOOD 

1 

1  ) 

9/ 

3/64 

0102 

07 

30 

21 

4b 

5N 

1  58 

1  0 

42W 

803. 

9 

00 

GOOD 

1 

20 

V/ 

J/o4 

0200 

46 

00 

21 

48 

6N 

1  bb 

1  0 

4bW 

843. 

4 

50 

GOOD 

1 

21 

9/ 

3/6 '« 

0J00 

50 

00 

21 

48 

8N 

1  58 

1  0 

4  9W 

850. 

9 

00 

GOOD 

1 

?d 

9/ 

3/6  4 

0904 

1  1 

40 

21 

47 

5  7N 

1  5B 

1 1 

4  W 

B8I  . 

9 

OO 

GOOD 

1 

23 

9/ 

3/64 

1  1  00 

4  4 

70 

21 

47 

4BN 

158 

1 1 

2bW 

854. 

4 

50 

CAPS 

1 

24 

9/ 

3/64 

1  304 

4  2 

70 

21 

47 

5riN 

1  58 

1 1 

low 

863. 

4 

50 

GOOD 

1 

25 

9/ 

3/64 

1500 

40 

70 

.21 

47 

35N 

1  58 

1 1 

23W 

854. 

4 

50 

GOOD 

1 

26 

9/ 

3/64 

1  700 

29 

39 

21 

48 

2N 

158 

1  0 

4  4W 

845. 

4 

50 

GOOD 

1 

27 

9/ 

3/64 

1800 

30 

40 

21 

47 

ON 

158 

10 

4  1  W 

858. 

4 

50 

GOOD 

1 

2H 

9/ 

J/64 

1900 

4  7 

20 

21 

48 

7N 

1  38 

10 

35W 

820. 

7 

73 

GOOD 

1 

29 

9/ 

3/64 

2000 

46 

29 

21 

48 

YN 

158 

10 

3  7  W 

761  • 

4 

50 

GOOD 

1 

30 

9/ 

3/64 

2101 

1  0 

oo 

21 

47 

ON 

1  58 

1  0 

54  W 

730. 

1  3 

bO 

GOOD 

1 

31 

9/ 

3/64 

2200 

38 

60 

21 

47 

59  N 

1  38 

1  0 

50W 

860. 

4 

bO 

GOOD 

1 

32 

9/ 

3/64 

2300 

38 

40 

21 

4b 

3N 

1  58 

1  0 

4  5W 

860. 

1  1 

2b 

GOOD 

1 

33 

9/ 

4/64 

0  300 

26 

60 

21 

47 

ON 

1  38 

10 

4  6W 

802, 

4 

50 

GOOD 

1 

34 

9/ 

4/64 

0100 

25 

40 

21 

48 

1  N 

1  38 

1  0 

57W 

774. 

1  1 

23 

GOOD 

1 

35 

9/ 

4/64 

0200 

22 

45 

2 1 

48 

YN 

1  50 

1  0 

43W 

802. 

4 

50 

GOOD 

1 

36 

9/ 

4/64 

OJ04 

20 

50 

21 

4b 

1  5N 

158 

1  0 

39W 

788. 

9 

00 

GOOD 

1 

37 

9/ 

4/64 

0400 

4  1 

20 

2  1 

48 

1  6N 

138 

10 

59W 

766. 

13 

30 

GOOD 

1 

30 

9/ 

4/64 

0  500 

29 

80 

21 

4b 

2N 

1  38 

1  0 

53W 

817, 

4 

30 

GOOD 

1 

3 '7 

9/ 

4/64 

0604 

53 

40 

21 

47 

22N 

1  58 

1  3 

6W 

831  . 

9 

00 

GOOD 

1 

40 

9/ 

4/64 

0702 

03 

20 

21 

47 

53N 

1  38 

10 

36  W 

729. 

4 

50 

GOOD 

1 

4  1 

9/ 

<■  '64 

ObO  1 

01 

20 

21 

4b 

1  N 

1  38 

1  0 

24w 

883, 

9 

OO 

GOOD 

1 

42 

9/ 

4/64 

0900 

00 

60 

21 

4b 

l  OH 

1  38 

1  0 

34W 

84  1  . 

4 

30 

GOOD 

1 

43 

9/ 

4/64 

1  000 

34 

00 

21 

4b 

4N 

1  38 

1  0 

42W 

773. 

1  3 

50 

GOOD 

1 

44 

9/ 

4/64 

1  7C  1 

1  4 

40 

21 

47 

44N 

1  30 

10 

4  6W 

801  . 

9 

00 

GOOD 

1 

45 

9/ 

4/64 

1  HO  1 

49 

60 

21 

48 

0N 

150 

9 

50W 

831  . 

4 

50 

GOOD 

1 

46 

9/ 

4/64 

1  ^>00 

38 

60 

.  \ 

4b 

4  N 

1  58 

1  0 

26W 

T04. 

4 

50 

GOOD 

1 

47 

9/ 

4/64 

2000 

37 

40 

21 

4b 

1  ON 

1  58 

1  0 

1  6W~ 

—  589. 

4 

50 

GOOD 

1 

4« 

9/ 

4/64 

210  1 

36 

00 

21 

47 

37N 

1  38 

10 

29W 

5d6. 

9 

00 

GOOD 

1 

4  0 

9/ 

4/64 

2200 

1  9 

ao 

21 

4b 

6N 

1  38 

9 

56W 

831  . 

1  3 

50 

GOOD 

1 

50 

9/ 

5/64 

0103 

57 

20 

21 

47 

57N 

1  50 

10 

1  9W 

766. 

9 

00 

GOOD 

1 

51 

9/ 

6/64 

0300 

49 

30 

21 

47 

2bN 

1  58 

10 

29W 

784. 

9 

00 

GOOD 

1 

52 

9/ 

5/64 

0400 

43 

1  0 

21 

47 

3bN 

1  38 

10 

1  3W 

659. 

9 

00 

GOOD 

1 

5  J 

9/ 

5/64 

0502 

46 

80 

21 

47 

49N 

1  38 

10 

ow 

624. 

9 

00 

GOOD 

i 

54 

9/ 

6/64 

0600 

2  b 

60 

21 

47 

49N 

1  38 

1  0 

26W 

663. 

22 

00 

GOOD 

l 

55 

9/ 

5/64 

0700 

23 

65 

21 

47 

3  l  N 

158 

1  0 

3W 

591  . 

10 

00 

GOOD 

1 

56 

9/ 

5/64 

1200 

32 

00 

21 

47 

43N 

1  50 

10 

43W 

667. 

9 

00 

GOOD 

1 

57 

9/ 

5/64 

1  300 

31 

20 

2i 

47 

55N 

1  38 

1 1 

36W 

604. 

10 

00 

GOOD 

1 

5b 

9/ 

5/64 

1402 

12 

.0 

21 

37 

52N 

150 

0 

57W 

671. 

18 

00 

DUD 

1 

60 

9/ 

5/64 

1600 

48 

.70 

21 

47 

50N 

1  30 

10 

59W 

615. 

18 

00 

GOOD 

t 
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APPENDIX  lla 
GM  BUOY  LOG 


GM  OAHU  SHOTS 

The  following  is  a  log  taken  from  the  GM  DRL  buoy  roll  chart.  The  column 
headings  for  station,  date,  and  shot  number  are  self-explanatory,  as  are  those 
for  hours  and  minutes.  The  fraction  of  a  minute  and  amplitudes  are  all  measured 
in  hundredths  of  an  inch.  The  last  number  is  a  calibration,  the  number  of 
hundredths  of  inches  per  minute.  Thus,  the  first  return  arrived  at  0845  and 
(z757  )  x  seconds,  lasted  for  ^  x  60  seconds,  and  had  an  amplitude 

of  1.75  inches. 


ALEUTIAN  TI  SHOTS 

A  reproduction  of  the  receiver  log  is  shown.  In  contradistinction  to  the  Oahu 
shot  log,  the  column  labeled  TIME  indicates  the  nearest  time  fiduciary  mark 
preceding  the  onset  of  the  arrival  signal.  This  time  is  given  in  hours,  minutes 
and  millimeters.  The  column  labeled  DUR  gives  the  duration  time  of  the  arrival 
signal  in  millimeters  and  SOFARR  gives  the  time  elapsed  in  hundredths  of 
millimeters  from  the  fiduciary  mark  to  the  arrival  of  the  axial  ray.  The  con¬ 
version  factor  is  60. 5  millimeters  per  60  seconds.  Thus,  for  shot  23  the  axial 
ray  arrives  at  2405  hours  plus  (69/60.  5)  (03  +  30. 00)  seconds,  and  the  duration 
is  (60/60.  5)  (73)  seconds. 


ARRIVAL  TIMES  OF  TI  ALEUTIAN  SHOTS  AT  GM/DRL  BUOY 


Stn 

Date 

Shot 

Time 

Dur 

SOFARR 

No 

mo  dy  yr 

no 

hr  mn  mm 

mm 

mm 

11 

090264 

23 

240503 

73 

3000 

11 

090464 

22 

243154 

76 

7150 

11 

090764 

11 

062113 

76 

3200 

11 

090764 

13 

191700 

74 

1900 

11 

090764 

26 

235730 

64 

4650 

11 

090764 

25 

032400 

64 

1600 

11 

090864 

24 

065034 

72 

5250 
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AWW1VAL  TIMES  GM/DRL  OAHU  SHOTS  AT  G~/DHt-  BUOT 


4  TN  OATE  SHOT 

NO  MOnVYWNO 
I  1007090  1640101 
I  1 00  7000 1640201 

I  10070  <01640501 
I 10070901640601 
I 10070901640701 

II  Of.  0901640801 
I  IOC  '902640901 
1100'.  ^02641001 
I 1007090264 1201 
I  10070 902641 301 
1  10070*026-*  1501 
1 10070902641601 

1  1 00  7090364 1 90 1 
I  100  70903642001 
1  10(J709'iJ642l01 
1 10070903642201 
1 10070903642401 
1 10070903642501 
1 10070903642601 
1 1 007090364270 1 
1 100 70903642301 
1 10070903642901 
1  1 00  70 90 364 300 1 
1 10070903643101 
1 10070903643201 
1 10070904643301 
1 10070 904643401 
1 1 007090464330 1 
1 1007090464360 1 
1 1 007090464370 1 
1 10070904643801 
1 1007090464390 1 
1 1 007090464400 1 
1 10070904644 101 
1 1007090464420 1 
1 10070904644301 
1 10070904644491 
1 10070904644561 
1 1 0070904644*01 
1 10070904644701 
1 10070904644801 
1 10070904644901 
1 10070905645001 
1 10070905645101 
1  10070905645201 
1  100  70905645301 
1  100  70905645401 
1  1 007090554550 1 
I  10070905645601 
1 10070905645701 
1 10070905646001 


TIME  AMPQUW  IN/ 
HRMN  IN  IN  |N  M|N 
0*145060  175032  237 
I044IH5  143025  236 
1645054  12  1029  230 
1844238  173030  239 
204/154  138031  237 
2245034  177033  236 

0047034  183034  239 
03440-J6  179026  2JH 
1 1 45029  1 78028  23d 
134  /165  165031  238 
1 744 193  1 35025  236 
234 Z034  172030  23/ 
0146085  173036  238 
0244237  172027  238 
0.344253  174034  238 
094809H  110023  238 
1  148229  1 78027  236 
1 34 4 2 1 H  178033  238 
1744430  155026  605 
1844179  172030  238 
1944243  1 /30 J6  238 
2044238  179032  238 
2145095  183034  238 
2244532  168026  605 
2344210  190040  238 
0044162  178035  238 
0144160  182036  238 
0244148  166030  238 
0348141  172033  238 
0444234  174034  239 
0544183  150025  239 
0649040  160032  239 
0746071  165029  238 
0845061  180029  238 
0944060  165028  236 
1044194  175030  238 
1  7451  I  0  0620  15  2 Jo 
1845250  135028  236 
19442QB  1 76030  238 
2044205  156030  239 
2145198  166030  239 
2244134  175034  239 
0 1 48050  184034  239 
0344255  173023  239 
0444230  183026  239 
0546244  IV20J2  23V 
0644157  132021  23 V 
0744150  100025  239 
1244186  098025  23V 
1344205  060019  23V 
1644294  178034  23V 
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APPENDIX  Mb 

MILS  STATION  RECEIVED  DATA  LOG 


The  following  listing  shows  the  recorded  date  at  MILS  stations  Eniwetok  and 
Midway  of  the  GM  DRL  Oahu  shots. 


The  column  headings  are  self-explanatory. 


ARRIVAL  TIMES  GM/UWL  OAHU  SHOTS  AT  MIDWAY 


STN 

NO 


datf 

SHOT 

time 

DUR  S  ^AR  AHR1 

VAL 

SEC 

SHOT 

M30YYR 

NO 

HPMN 

SEC  CHI 

CH2 

CH3 

CH4 

BRG 

9 

264 

1  8 

12325 

1274 

377 

9 

36* 

19 

10124 

1403 

506 

9 

364 

20 

1  0223 

1  188 

290 

9 

364 

21 

1  0323 

1  228 

331 

9 

364 

22 

1  0926 

144  1 

543 

y 

364 

24 

11327 

1  1  53 

255 

y 

364 

26 

1  1  523 

1  133 

237 

y 

364 

26 

1  1  723 

1021 

124 

9 

364 

28 

l  1923 

1  20  1 

304 

V 

364 

29 

12023 

1  188 

291 

9 

364 

30 

12123 

1428 

530 

9 

■>64 

31 

12223 

1112 

216 

9 

364 

32 

12323 

1113 

216 

9 

464 

46 

1  1923 

1116 

217 

o 

464 

47 

12023 

1  1  04 

207 

9 

464 

46 

12124 

1097 

1  97 

9 

464 

49 

1  2223 

938 

041 

9 

364 

-30 

10126 

1  296 

398 

9 

564 

51 

1  0323 

1  229 

331 

9 

564 

66 

1  1  223 

1046 

151 

9 

664 

67 

1  1  323 

102  7 

1  28 

9 

664 

60 

1  1623 

1215 

318 

i\ 

264 

18 

12351 

420 

4  76 

342 

044 

9 

•»64 

20 

1  024R 

952 

1007 

873 

555 

9 

364 

21 

1  0348 

992 

1  050 

91  J 

596 

9 

364 

2? 

1  0952 

610 

669 

536 

215 

9 

364 

24 

1  1  352 

91  4 

971 

838 

51  7 

9 

364 

25 

1  1  548 

899 

957 

820 

503 

9 

364 

26 

1  1  748 

788 

844 

706 

391 

V 

364 

28 

1  1948 

965 

1023 

888 

572 

9 

364 

29 

12048 

955 

1012 

879 

556 

9 

364 

32 

12348 

876 

935 

798 

481 

9 

464 

47 

12048 

869 

924 

792 

475 

9 

464 

46 

12149 

868 

926 

783 

4  75 

9 

464 

49 

1  2248 

704 

r  58 

624 

304 

V 

S64 

51 

10348 

988 

1  046 

910 

591 

9 

564 

52 

10448 

926 

984 

849 

532 

9 

564 

53 

1  0550 

969 

1  026 

892 

571 

9 

564 

5* 

1  0649 

162 

218 

086 

9 

564 

55 

1  0748 

735 

791 

659 

339 

9 

564 

57 

1  1  348 

788 

846 

710 

393 

9 

564 

60 

1  1648 

984 

1040 

907 

387 

PIS  UG 

Ot)  06V  A  GOAO 
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APPENDIX  III 

SOFAR  VELOCITY  PROGRAM 


The  following  is  a  listing  of  the  program  used  to  convert  source  and  receiver 
logs  into  average  sound  velocities.  The  program  mechanics  are  briefly  described 
in  Section  €.  The  coding  is  done  in  IBM  7040  Fortran  IV,  and  the  program  was 
run  on  an  IBM  7040  computer. 


C  PROG-'AM  T U  COMPUT-  AVERAGE  SOUND  VELOCITY  FROM  SOURCE  AND  STATION  LOGS 
;  .•■'rtrNNAN  1/16/65 

DIMENSION  1MOHOO  )  , SDAV!  IOOI  .SHW!  I  00  )  •  SMN  !  I  OC  )  .SSEC<  tOO  I  «  SLAT  I  1001 

1  . SLNG! 1 OO ) .SOEPTH! 1 00  I . SWT  ! 100)  .SPEHI I  00) .RSCI A  > .Nl  I  10) .N2I60) 

2  .RCLATIA ) ,RCLNG!A ) . ROEPTH ( A ) . L AOEL < 56 ) 

=^AL  N 

Rc  ao  2.N! .N? 

2  r OPVAr < 1 OA 1 , 35A2/25A2) 

Rr.Ao  3.N.S.F.W.G00D 

3  FORMAT  t  A  A  1  • A A  ) 

REAP  Z''  «LA0FL 
FORMAT  (  1  ->A6,A2  ) 

A  RE  AO  5.  I  SRC.  1  SHOT. |M0I  I  SHOT ) .sUAYl | SHOT ) .SHWI  | SHOT ) .SMN I  I  SHOT ) . 

1 SSEC < ISHOT ),SLATOG.SLATMN.SLATSC.SLATNS.SLNGOG.SLNGMN>SLNGSC. 

2  SLNGE'V.SOEDTH)  I  SHOT  )  .SWT  !  1  SHOT  I.SPERI  ISHOT  ) 

5  FORMAT |  |? <2 | 3. >  X .F2«0. AX.2F2.0.F6.2. 1 X. 3F3. 0 . A  1 .2X.3F3.0.AI  ,3X. 

1  FA.O.FT, ?, l x. AA 5 

|P( ISHOT. EQ.100IG0  TO  10 

SLATRDs  .  I  7  1532925 E-  1  *  I  SLA  TOG  +SLA  TM«>:/60.  +SLATSC/3&00. > 

SLAT ( ISHOT )=SLATR0~ .337267 I 2t-2»S INI  2. *SL A TRO > ♦ »56B7E-b*S IN  I A.* 

I  SLATRO) 

i  f  <  slatns. eq  *  s ) slat i ishot)»-slat4 ishoti 

SLNG  ! I SH0T )  =  « 1 7AS32V25E- I  *  I SLNGDG  ♦SLNGMN/60*  ♦SLNGSC/3600. > 
IF(SLNGEW.EQ.V)SLNGI ISH0T)«  -SLNG! ISHOT) 

GO  TC  A 

10  no  u  i=i,2 

READ  1  1  . PLAT , RCLNG (  1  ) .ROEPTH 1 1  ) 

11  FORMAT  (?F15«r'»F10.0) 

RLAT*RLAT  -,337267|2E-2»S|NI2,*KLAT)  +.5687E-5*S|N! A.ARLAT ) 

RCLAT1  |  ) *RL A  T 
RCLAT( | *2 ) «RL AT 
RCLNG  <  1+2  > «RCLNG  < I  ) 

1A  ROEPTH! | F2 ) «ROEPTH ( | ) 

M*0 

15  READ  12 . ISTN.ROAY. 1  SHOT . JSWC. I  SETUP. ICHAN.MHR.RMN.RSCST.DUHMIN. 

I  OURSE C.RSC.  1  AMP.  DQR  I  SE  ,  DtJtJKGD .  A T TEN ,  GZERO  .  AT2F.H0 

12  FORMAT! l2.bX.F2tO.2X. 14. 12 , I* . 1 2.SF2.0, QFS.2.4X. I 2.2F3. 1 »F2. 0.F3.0 

1  .F2.0) 
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i  i  ^ '  n  •  g  t  « 2  o ) go  ro  10 

•••  -  '*  U 


l-'  issc''T.c,r,i?,jf;  <  ichani  ) Fuuot  =  t>0 • 

1  I  \’r  =  flf>4  "> ')»  »  < -VI.) A  Y  -  SDA  V  <  I  SHOT  )  I  +3600,  *  1  HKW-SHP  (  I  SHOT  )  ) 

1  jvn-'j^NI  If.HO'  )  )  FPbCI  ICMANI-SSEC  «  I  SHOT  )  ♦ruOGfc 

;'l4t  =prLar  <  i  i-M6N  i 
xLN'is^'I.NJ  <  I  CHAN  ) 

-All  go  t  -'la  t  .wlng.slati  i  jhj  r  i «  slng  i  i  shot  >  «  mange  «e*.eni 
•’•vlats  plat  +  slat < 1  shot n 

H\  1  lu^-.V!  OOG.  *  <  .9^63200 4  7  l6aJ49AE-2*coS<2»*AVLAT  ) - •  3S49E-3* 

I  ."><5  (4,»  A VI. AT  »  ) 

■;.*-\Gr  i  u~-,*pange 

C -*&•,'  =  PAW  /T  |  ME 
^<v:^A\r,r*  ,if1 1 
AT  =  P<-f  (  1  )-  QSC  <  2  ) 

!-'  IWSC  U  ).c3.0».  ANP.RSCI2  »  •tU.O.  »OT*PSC « 3 > -RSC  « 4 > 
tPL-0.)9  I  r,“*l'J>3\ND+ATTtN  -GZEH0-ATZEQ0F80* 

p«VsfAAV(  I  GhOT  )  ♦  1  , 

I v | N  =  S  VN  <  I ' H  3  T  )  +1. 
t  HP  =  'MD  <  I  SHOT  )  «  |  • 

:r  «*frn(v.40»,ra,i  i print  20 

?0  F  -'OVAT  IIHI  I 

lr  (VCD('«,«PI.E0.I  IPPINT  22.LA0EL 

P.9  IN  T  24 . JSPC* IoTN,  I  SHOT  ,  |MO<  I  SHOT  I  .  Nc  «  I  DAY  )  ,  N2  «  I  HR  )  .  N2  C  I M I N  )  * 
1  OUAK.SDEPTHI  I  SHOT  I  , ROEP  Tm (  I  OMAN)  .DT.SPL 

PUNCH  24, JSHCi ISTN, ISMOT . I  Mg ( I  SHOT > ,  N2  <  I  DAY  >  «N2  <  I  HR  )  •  N2  <  I  M 1 N  >  • 

|H<y,C0AR,SDEPTM( 1SH0T|,R0EPTK< ICMAN) tDT.SOU 
24  F0HA’AT<  4  1  3,  |H/,  A2,4H/6«  ,  2  ',2  ,  F  V  .  2  ,  F9  •  2  «  2F6»  0  ,  F  7. 2  «  F6»  0  *  1 43  ) 

GO  TO  |« 
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c 
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I 
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APPENDIX  IV 

AVERAGE  SCIFAR  VELOCITY  LISTINGS 


The  following  are  listings  of  the  output  from  the  sound  velocity  program 
described  in  Section  8  and  Appendix  m.  Average  velocities  are  given  for  the 
GM  DRL  deployable  monitoring  buoy  station  for  both  the  Oahu  and  T1  2376-  lb 
shots.  The  column  headings  are  self-explanatory  (ZSOURCE  and  ZREC  are 
depths  in  meters. 

AVERAGE  SOUND  VELOCITIES  TI  ALEUTIAN  SHOTS  IN  GM/DRL  BUOY 


SCR 

RC 

Shot 

Date 

Source 

Range 

C  Avg 

Z Source 

Z  Rec 

No 

No 

No 

mo  dy  yr 

Time 

km 

m/sec 

m 

:a 

2 

11 

23 

090264 

2311 

4805.93 

1466.70 

91. 

732. 

2 

11 

22 

090464 

2336 

4863.44 

1468.54 

91. 

732. 

2 

11 

11 

090764 

0527 

4803.48 

1468.06 

91. 

732. 

2 

11 

13 

090764 

1824 

4696.37 

1468.08 

91. 

732. 

2 

11 

26 

090764 

3305 

4649.46 

1468.49 

91. 

732. 

2 

11 

25 

090864 

0230 

4693.80 

1468.53 

91. 

732. 

2 

11 

24 

090864 

0557 

4746. 58 

1468.47 

91. 

732. 
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AVERAGE  SOUNi)  VtLOC 1 t I E S  G'l/OPL  OAHU  SHOTS  TO  G*/DWL  MOOT 


SBC 

BC 

SHOT  DATE  source 

Wi.MGl 

v-  *W5  ZSOUUCE 

2PEC 

NO 

NO 

NO 

N  ,)V  YK' 

T  1  «£ 

M 

M 

1 

1  1 

1 

'/<">1  /O* 

OaO  1 

»2-i«0* 

1  *  O 1  .  39 

•*20. 

fja* 

1 

1  1 

a 

*/01  /6« 

1000 

392  >.  73 

1  «  7  7.  M2 

72  0. 

732. 

1 

1  1 

n 

9/0  1  /b« 

1600 

3930. 3  £ 

1 *60.3* 

6*9. 

732. 

1 

1  1 

t> 

J/P 1 /6A 

-1U00 

3927.31 

1  *on,u9 

732. 

>32. 

1 

<  1 

7 

J/Al /6A 

200  3 

3926. *9 

1*83.03 

707. 

732. 

1 

1  1 

H 

'/O 1 /6A 

c.20  1 

3->26.  77 

1*80. 1 7 

676. 

732. 

1 

1  1 

9 

9/02/6* 

0002 

3U2 7.09 

1 «  79,63 

9*7. 

732. 

1 

1  1 

1  0 

9/02/6* 

0500 

3926.01 

1*  T9.21 

711. 

732. 

1 

1  1 

12 

9/02/6* 

1  100 

3926. 73 

1  *  79.  77 

8  72. 

7 32. 

1 

1  1 

1  3 

9/02/0* 

1303 

3*26.3 1 

1479.61 

601  . 

32. 

1 

1  1 

1  5 

9/02/6* 

1  700 

3 -*£3.00 

I*  7  ,.,41 

771  . 

Tj£  • 

1 

1  l 

IB 

9/02/6* 

2  302 

3925.  -96 

1  *  79.  31 

750. 

732* 

1 

1  1 

1  V 

9/03/6* 

0102 

3  92  3.96 

1  *  ft* 33 

603  • 

73c. 

1 

1  1 

20 

9/03/6* 

0200 

J92r>»  1  0 

l *  79.46 

8*  -• 

’32. 

l 

1  1 

21 

9/03/6* 

0300 

3926. 10 

l *  / 9,4* 

850. 

’32. 

1 

1  1 

22 

9/03/6* 

090* 

3926.6* 

1  *  79. 90 

861  • 

732. 

1 

1  1 

2A 

9/03/6* 

1  30* 

3926. 78 

1*  78,99 

863. 

732. 

1 

1  1 

25 

9/03/6* 

1500 

3927. *9 

1  *  79.69 

85*. 

732. 

1 

1  1 

26 

9/03/6* 

1  700 

J926.Pb 

1479.66 

8*5. 

732, 

1 

1  1 

27 

9/03/6* 

1H00 

3926.68 

1*  79.20 

658. 

732. 

1 

1  1 

2d 

9/03/6* 

1900 

3925. 76 

1*79,13 

820. 

732. 

1 

1  1 

2<? 

9/03/6* 

2000 

3923.81 

1*79.32 

761  . 

732. 

l 

1  1 

30 

9/03/6* 

2101 

392  7,2  1 

1  *  79,  76 

730. 

732. 

l 

1  1 

31 

9/03/6* 

2200 

392»>»2  3 

1*79.28 

860. 

732. 

1 

1  1 

32 

9/03/6* 

2  300 

J9c  6.07 

1*  79,00 

160. 

73c. 

l 

J  1 

33 

9/0*/6« 

0000 

392  7.01 

1*  79.62 

802. 

732. 

1 

1  l 

3A 

9/0*/6« 

0100 

392  6  »  *  0 

1*  7d,9| 

77*. 

732. 

1 

l  1 

36 

9/04/6* 

0200 

3923.96 

1*7«.  73 

802. 

732. 

l 

1  1 

36 

9/04/6* 

030* 

3925.7* 

1  *  76.60 

"788  . 

732. 

1 

37 

9/0*/6* 

0*00 

3926.2* 

1*77. 39 

766. 

732. 

1 

l  1 

38 

9/04/6* 

0500 

3926.29 

1  *  76,  1  9 

81  7. 

732. 

l 

1  1 

39 

9/0*/6* 

060* 

3930.2* 

1 0  79,00 

831  • 

732. 

l 

1  l 

AO 

9/0*/6« 

0702 

3926.*9 

1  *  79.0  7 

729. 

732. 

l 

1  1 

A  1 

9/0*/6* 

0B01 

3923.37 

1  *  79.01 

883. 

732. 

l 

1  1 

A2 

9/0*/6* 

0900 

3923.69 

1  *  76.67 

8*1  • 

1 32  • 

1 

l  1 

A3 

9/04/6* 

1003 

3923.98 

1*79,21 

7/3. 

732. 

1 

1  1 

A  A 

9/0*/6* 

1  701 

3926.37 

1  *  79,  7V 

HOI  . 

732. 

l 

1  l 

A5 

9/0*/6* 

1801 

3924.60 

1*79,07 

831  • 

732. 

1 

1  1 

A6 

9/0*/6* 

1900 

3923. 3  7 

l*  79,21 

70*. 

732. 

1 

1  1 

A7 

9/0*/6* 

2000 

3925.23 

1*76.93 

589. 

732. 

1 

1  1 

AB 

9/0*/6* 

2101 

3926.0* 

l*  79,46 

5B6. 

732. 

1 

1  1 

A9 

9/0*/6* 

2200 

392*. 78 

1*76.91 

83 1  • 

732. 

l 

1  l 

50 

9/05/6* 

0103 

3925.30 

1*76.33 

766. 

732. 

1 

1  1 

51 

9/05/6* 

0300 

3926. 1 7 

I  «  76»'9A 

78*. 

732. 

1 

1  1 

52 

9/05/6* 

0*00 

3923. S2 

!«  78.  78 

659. 

732. 

1 

1  1 

53 

9/05/6* 

0502 

3923. 1 3 

1*  70.  70 

62*. 

732. 

1 

1  1 

3* 

9/05/6* 

0600 

3925.79 

1*79. 30 

663. 

732. 

1 

1  1 

55 

9/05/6* 

0700 

3925.52 

1*79.06 

591  . 

732. 

1 

l  1 

5ii 

9/05/6* 

1200 

3926.31 

1*79,01 

667. 

732. 

1 

1  1 

60 

9/05/6* 

1600 

3926.61 

1*78.91 

615. 

732. 
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APPENDIX  V 

SOURCE  LOCATION  PROGRAM 


The  following  is  a  lift  ting  of  the  GM  DHL  source  location  program  configured 
especially  for  the  University  of  Hawaii  ML2  data  format.  The  coding  is  an 
IBM  7040  Fortran  IV,  and  the  program  was  run  on  a  1040  with  a  real-time 
storage  clock.  The  program  is  described  in  greater  detail  in  Section  6. 


C  PROGRAM  TO  COMPUTE  SOURCE  LOCATION  or  UNDERWATER  ACOUSTIC  EVENTS 
C  6  V  Li  AST  SQUARES  CRITERIA 
C  **  vLLFNNAn  JAn  I96S 

DIMENSION  HLATI-J9  I.WLNGt  V9I.RCAVI  99  >  •  S I  NUTA  125  >  •  COiOT  A  1 25 »  •  T  ABl  25 » 

1  .SLAT  I.SLNGI2S  I. SC AVI  2SI.JI2V.  .CRAOI2SI*  I  DAY  1261  .  1*0(251 

2  •  NSTN I  ?*5 1  »  PANGS  <251 

External  ipcl* 

p-al  n 

30^T<v>  =  v.,337267|2E-2*SIN<2.*V)  •  ,56B7E-5*S INC  a. *y » 

G3CN(''|  =y-,33726712E-2*SIN(2.»Y»4.56B7E-5«S!NIA.*YI 
IT|«E«t>=4096*MOC< INT(T/60. »  *  I  0  >  4262IA4*  *001  INTI T/ftOO. > .61 
I  -H6777?|6*MOO<INT IT/3600.  I.IOI  410737*IB2A*  INT I T/ 36000. » 

pal  IV >=6378380.* I .9983200*74. I 6B3494E-2 *COS (2 • • Y ) - • 3S49E-S*C0S ( 4 . • 

1  v  >  » 

•JEAo  2.N.S.E.V 

?  FORMAT | 4A1 » 

C  INPUT  STATION  OATA 

A  -(HAD  5. I STN.RLATI 1STN) .RLNGI 1 STN) .RCAVI ISTN) 

B  FORMAT  I  I2.FI3.B.F15.8.FI0.2) 

Ic( ISTn,c0.Q91G0  TO  10 

RlATUSTN»  =  SIGNCGOCNIABSCRLAT  1 1  STNI  I  )  .PLAT  I  ISTNI  I 

C.1  TO  A 

c  Input  trial  parameters 

10  PEAO  I ?. V2EPO.X2EPO 
12  F0pmat«?F|0.2» 

X2ER0=.0I 74S33*XZER0 
Y2EPO* . 0 1 7*b33*v2EPO 
GOTI*E=I ABSl ibclxi 

C  INPUT  STATION  ARRIVAL  TINES 
I  *0 

I*  Ul*l 

17  FOR**AT  <  I2.6X. I2.5X. I2.4X. 12. I J.2X. I2.3X.F3.0.2X.F5.3.9X. I2> 

READ  I  7. JShOT . I  DAY  I  I  ) . INOII  ». |YR. I  HR. I *N . T SC. TF ACT . ISTN 
IF» JSMOT.NE. ISHOT »GO  TO  20 
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I  H  Jl'Tin  xC*.  AT  (  ISTNI 
SENG  I  |  )  rPL-jr,(  IST'J  ) 

C A V  I  1  )=PCAVI  ISTNI 
N'TNI  I  |3TN 

'WII:  3600.  *FLOAT<  |mmi*60.*FlUAT»  I  NN  )  ♦  T  SC  *  TF  AC  T*  t  ,£>A00.  »TLOAT  ( 

I  | ^AV (  I  )  -  | DA  V  < I  |  | 

I ^HOTrj^HOT 
03  TO  t* 

^  i  too  to  i  a 

>cxrt'^o 

v  I  v  *  -  p-> 

f»T«(|  1-1699. 

3  S’’ APT  Irr-PATJOM 

CO'JO'jTC  C3t  c  =■  I  C  I^NTS 
1  T\"’  =  3 

’A  |TVI:|tn(1  *  | 

’I  1  =  1  «  |NAX 

viLL  GCOISlAT  «  I  I .SLNG«  1  »  *  V .X. RANGE (!  I.SI N8TAI  |  ) .COS8TAC I  )  ) 
Avl/iTj.N.i  v*SLAT  mi 
CPAn« I )xSCAVt 1 l/OAPIAVLAT J 
•»0  ?lll«CP«f)ll)»ITARt||-T|  -RANGEtl) 

-ALL  LSTSQJISINHTA.COStJTA.LRAO.D.OX.JY.DT  iME)<  1  NAX  i 

Xx  X-PX 

YxY-PY 

TsT*DT 

yxAbsiavodii  ,b70  7963-v.6.2H3lt'53»  ) 

|E|v.ST.-».|A|S927)Yx6.283I8b3-Y 

Yx| 

X  =  4V"(X  BXT) 

lcl»«S(x  i.GT.3. 1AI 9927 ix=SIGN 16. 2831 8brt-Aascx>  .-X  ) 

I c  (  I  TNO.GT  .  .?•>  1 00  TO  32 

|«t»ars|OXI.C,T., 158^-5. OR.  ABS«DY).GT.I58E-b>G0  TO  2b 

C  OUTPUT  PENULT*; 

15  PP  (NT  79,|«>0T 

3b  FORMAT  «  I  Hi  X  AJMV.EAST  SQUARES  pOUPCt  LOCATION  OF  T.I.iMOT.lJ  // 

1  2 7m  APPIVAL  tines  at  stations  /a oh  station  date  gct 

2  RESIDUAL  /39H  NUN0EO  NO  0 Y  YH  HRNN  SEC  NETENS  /  I 

ro  -XH  I  X  I  ,  |  VAX 

P=-SPuTs*37inO0»*OI  I  I 

jT|»3t|T  I xr (TAPI  I  I  I 

Str  r  £V6n|TARI I ) .60. I 

10  POINT  lp> •  NSTN  (  I  )  »  I  NO (  I  I  •  I  DAY  (  I  I  •  |  YR »  JT  I NE •  SEC •  RESOUT 
36  FORMAT 1216*21 IM/. 12) .2X.A4.F6. 2.  F8.0) 

PLATNSiN 

IE| Y.LT.O. |PLATNS»S 
OV.ATxS7,298779S»G00T  lAOSIT  »  » 

ILATDGxOLA? 

II.ATMN*  60  «  ®  AMOO  C  PL  AT  «  I  .  ) 


\V 
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ll  AT'-Ct  AMODI 3600. *PLA t  .60.  > 
nj  Mr.c*  »i-r- 

I  "  (  v.L  r  ,  o.  )DLNGFW=W 

3|  Mr=  ?t:77'»R#  ATS  (  X  ) 

II  Mri'A-OLVA 

1 1  vr-VMKSO.  #  A  MOO  <  PL  NO.  .  1  •  ) 
lLMr-=.Cr  AV0'1<  "»600.*PLNG.o0.  > 
JT | vr- | T |M- (T ) 

Giro  -  AMOO(T.60.  ) 
lr'*S> 

RrS0UT  =6371000.*RES 


Point  AO 


40  FORMAT < 

1  NO. 

/ / / 6 1 H  SHOT 

/ 

DATE  GCT 

LAT 1TUDL 

LONG  1 TUDE 

RMS 

2 

64H 

MO  DY  YR  HRMN  SLC 

L>G  MN  SC 

DEG  MN  SC 

RES 

A  ITERATIONS  ) 

PR  I  NT  4  1.;  SHOT  .  1  MO  <1  >  .  IOAY ( 1  > .  I YR. JT 1  ML.  SEC.  IL  A  TOG.  ILATMNI.  1LATSC. 

1  PL  A  TNS.  1LNGDG.  1  LNGMN  .  1  LNGSC  .  PLNGLW .  RESOUT  .  I  TNO 
A  1  FORMAT  <2H  2 , 2  I  3 . 2  (  1  H/  •  1  2  >  ♦  1X.A4.F6.2.  1  4 . 2  I  3  «  A  1  .  1  6  '  2  1  3  •  A  1  .  F6 . 0  «  16) 
PUNCH  4  2.1  SHOT. I  MO  <  1  )  . 10AY  < 1  ) . ! YR . JT 1  ME • SEC . 1LATOG.  1 L ATMN . L L A TSC . 

1  PL  ATMS. ILNGOG. I LNGMN . 1 LNGSC . PLNGEW . BESCUT 
4  2  "OPMAT (2H  2 «  2 1  i »  2 ! 1 M/ « I  2  >  »  1 X  •  A  4  «  F  6 • 2  « I  4  -  2  I  3  •  A  I  «  I  5 «  2 1 3 «  A  1  .F6.0.26X  ) 
TLAP-ALr  (FLOAT  <  I  ABS  <  I  BCLK  )  )  -GOT  1  ME  >060. 

GOT  1\T:  1  ARS  (  IRCLK  ) 

PP.NT  4G.TLApSE 

44  F0PMAT(//1HH  COMPUTATION  TIME*  .F7.3.  4H  SEC  ) 

C  SET  UP  FOR  NEXT  SET 

I  DAY (1)= 1  PAY  <  1MAX+1  ) 

1M0<  1  >= 1 MC< I  MAX* I  ) 

I  *  1 

GO  TO  i 
ENO 
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APPENDIX  VI 

SOURCE  LOCATIONS  OF  248-lb  SHOTS 


The  following  Is  a  tabulation  of  probable  source  locations  for  the  ?I  248-lb 
shots.  Three  sets  of  figures  are  given  for  each  shot. 

a)  The  output  of  the  GM  DRL  program 

b)  The  radar-derived  position  of  the  plane  at  the  time  of  drop 

c)  The  output  of  the  University  of  Hawaii  epicenter  location  program 

Lines  a)  and  c)  differ  only  in  that  the  University  of  Hawaii  program  used  only 
the  MILS  data,  while  the  GM  program  used  additional  data  from  two  West 
Coast  stations. 

In  the  column  marked  Res-RMS,  the  GM  station  residual  is  expressed  as  rms 
meters,  whereas  the  University  of  Hawaii  residual  is  apparently  expressed  as 
the  mean  of  the  absolute  value  of  the  res' '  al  in  seconds.  There  is  no  residual 
expressed  for  the  radar  fixes.  Since  thebe  source  cards  are  in  the  same  format 
as  the  TI  source  logs,  the  extra  three  numbers  refer  respectively  to  shot  depth 
(meters),  charge  weight  (pounds),  and  local  water  depth  (meters). 
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SOUPCE  PO-,.1, 0N5  OF  T.,.  248-L8  ALEUT  1  AN  SHOTS  AS  COMPUTE  BY 


SHOT 

2  6b 
2  6b 

MO  DY  YP  HPMN  sec 
DATE  GCT 

9/17/64  2269  66.74 
9/17/64  2269  61.00 
2259  54 

DG  MN  SC 
LAT 1 TUUt 

51  27  22N 
51  35  00N 
51.5  N 

DtG  MN  SC  HES 
LUNGITUDt  WMS 

1 76  39  12w  2154. 
176  38  24»)  J8. 

1  '6.6  W  1,3 

2*8.  183 

2  66 
2  66 

9/17/64  231 d  1.62 
9/17/64  2314  69,00 
2314  54 

51  13  32N 

5 1  25  00N 
51 .3  N 

l'6  Je  1 3W  5061 • 

1 76  38  24W  38, 

176.6  W  3.5 

248.  1428 

2  67 
2  67 

9/17/64  2331  41,85 
9/17/64  2331  34,00 
2331  36 

Jl  7  27N 
51  15  00N 

51 .2  N 

176  37  52W  4268. 

1  '6  38  24W  J8. 

1 76.6  W  3,0 

248.  3480 

2  6B 

2  68 

9/18/64  -014-33,20 
9/17/64  2346  26.00 
2345  28 

60  Dti  29N 
51  05  OON 

61 ,0  N 

176  36  4 Vw  4057. 
176  38  2 4*  38. 

176.6  W  2,9 

2A8.  3850 

2  69 

2  69 

9/18/64  0030  27.16 
9/18/64  0030  37,00 

0  30  29 

60  59  3N 
60  36  OON 
51.0  N 

176  31  50W  7039. 

1 '6  38  24W  J8, 

176.5  w  2,7 

248.  3380 

2  70 

2  70 

9/18/64  0044  39,68 
9/18/64  0044  46,00 
0044  42 

50  «9  2N 
60  45  OON 
50,8  N 

1  '6  35  23W  6192* 

1  ^6  38  24r  J8, 

1  '6.6  W  2.2 

248.  4395 

2  71 

2  71 

9/18/64  0100  5.31 

9/18/64  0100  13.00 
0100  15 

60  4  8  UN 
60  J5  OON 
50.7  N 

176  35  Owl  1 137, 

176  38  24 W  38. 

1 76.6  W  2,9 

248.  6220 

2  72 

2  72 

9/18/64  0115  14,08 
9/18/64  0115  19.00 
0115  20 

5 J  28  39N 

50  22  51 N 
50.4  N 

176  36  45W  3926, 

1 76  40  OOW  38, 

176.6  W  2,6 

248,  6960 

2  73 

2  73 

9/20/64  2159  57,69 
V/20/64  2159  58,00 
2159  56 

61  53  22N 

51  52  27N 

51,9  N 

1 74  46  1 lw  2813. 

174  43  18W  38, 

174.8  W  1.1 

248.  156 
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SOU'  CE  POSITIONS  OF-  T  •  1  #  248-Lb  ALEUTIAN  SHOTS  AS  COMPUTtD  bY 


|.G”/nHL  PROGRAM  USING  U.OS  H,  DATA  PLUS  WEST  COAST  HYOROPHONtS 

2. NAVIGAT  IONAL  RADAR,  AS  REPORTED  IN  U.S.C4-6S  MEMO  DATED  l  OCT  I 964 

3. U.OE  HA*A||  PROGRAM  USING  MILS  HYDROPHONE  DATA 


MO  or  YR 
SHOT  OATt 
74  9/20/64 

74  9/20/64 

2  7b  9/20/64 
2  7b  9/20/64 


hrmn  sec 

GCT 

2100  14.99 
2  I  SO  lb* 00 

2iso  oa 

£236  JO • 32 
22J6  2r*00 
2236  33 


DG  mn  SC 
LAT I TUOE 
ol  3j  S9N 
Jl  32  SAN 
52.0  N 
b  I  4  1  2  JN 

61  43  I2N 
*J|.7  N 


DEG  MN  SC 
LONGITUDt 
I  7-j  26  3  7V 
I  7  b  23  24V 
I 7b*4  V 
I  7_ >  bb  394 
I  73  35  00V 
I7j.9  V 


RES 
RMS 
1 7 1  a » 
Jb* 

1*3 

49a  i  • 

3b. 

I  .3 


24B.  145 

24b.  tao 


2 

76 

9/20/64 

2132 

47.07 

51  36  28N 

176  14  | 7b 

31  34. 

2 

76 

9/20/64 

2132 

43.00 

51  40  3 ON 

176  15  24« 

38. 

240. 

165 

2132 

47 

bl.b  N 

1  '6.2  V 

1  .3 

2 

77 

9/20/64 

20*6 

33.66 

3 1  31  7N 

1 77  13  |4V 

1392. 

2 

77 

9/20/64 

2046 

46.00 

bl  40  JON 

1  77  04  OOW 

JB. 

24b. 

120 

2046 

57 

51  .5  N 

177.2  V 

1  •  1 

2 

7a 

9/20/64 

2032 

4  |  *b9 

bl  2b  46N 

177  2b  41V 

1070* 

2 

7a 

9/20/64 

2032 

30.00 

51  40  I2N 

1 77  24  24V 

Jb. 

24B. 

137 

2032 

44 

51  .4  N 

1 77,5  V 

l.» 

2 

79 

9/20/64 

2022 

34,  1  7 

bl  J9  31 N 

1 77  53  49V 

3469. 

2 

79 

9/20/64 

2022 

3b.  00 

b  l  Jb  OON 

1  77  47  4bV 

JB. 

248. 

147 

2023 

08 

51  .2  N 

1 77,9  V 

0.9 

2 

eo 

9/20/64 

1  94  3 

56.  7b 

bl  39  4 JN 

1 ’8  lb  17V 

8654. 

2 

ao 

9/20/64 

1  94  J 

59.00 

51  J3  42N 

1  7t_  lb  21  V 

38. 

248. 

210 

1944 

07 

51  .b  N 

1  76.3  V 

4.B 

2 

at 

9/20/64 

232  J 

1  9,  3b 

bl  29  b9N 

1  7b  46  1 bV 

5619. 

2 

81 

9/20/64 

2J23 

07,00 

bl  2b  42N 

1  7b  37  21V 

JB. 

248. 

320 

2323 

23 

51.4  N 

178,7  V 

2.7 

2 

62 

9/21/64 

-007-29.44 

bl  6  49N 

1  79  22  1  6V 

2571  • 

2 

82 

9/20/64 

2  352 

2  1.00 

51  12  24N 

179  16  00b 

JB. 

248. 

1815 

2362 

26 

51  .2  N 

179.4  V 

1  ,8 

2 

aj 

9/21/64 

00  lb 

2b.  JJ 

bO  36  JN 

1  79  3b  20E 

2461  . 

2 

83 

9/21/64 

ooia 

07,00 

51  12  OCN 

1  79  29  24E 

Jb, 

248. 

1095 

0018 

22 

51  .0  N 

1  7V.6  E 

I  .8 

2 

94 

9/21/64 

0040 

57,33 

bl  27~2 1 N 

1 7a  57  IOE 

625. 

2 

84 

9/21 /64 

0041 

02.00 

bl  19  b4N 

1 78  39  OOE 

38. 

248* 

210 

0040 

5b 

51  .5  N 

178.9  E 

0.5 

2 

16 

9/21/64 

0105 

19.69 

50  52  24N 

179  13  5BE 

3195. 

2 

16 

9/21/64 

010b 

1  0.00 

SO  59  3 ON 

179  10  36E 

38. 

248. 

3000 

0105 

20 

50.9  N 

179.3  E 

2.5 
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